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Abstract 


This brief review is devoted to the results of non-accelerator study of strong 
nuclear long-range interaction in the mass isotope effect by the modern method of 
low-temperature optical spectroscopy of condensed matter. This became possible after the 
discovery that the addition of one neutron causes global changes in the macroscopic 
characteristics of a solid. The observatian of an isotopic shift (0.103 eV) of the zero-phonon 
line of free excitons in the luminescence spectra of LiH (without strong interaction in the 
hydrogen nucleus) and LID (with strong interaction in the deuterium nucleus) crystals was 
the first and direct evidence of the long-range interaction of the Yukawa potential. indeed, in 
both crystals, the lithium ions, the proton and the electron are the same and, therefore, the 
gravitational, electromagnetic and weak interactions are the same, and the addition of a 
neutron, according to Yukawa, a strong interaction appears, the influence of which 
manifests itself in the isotopic shift. These experimental results demonstrate the 
neutron-electron binding energy (0.106 eV) which is in excellent agreement with the 
theoretical Breit estimate of 0.1067 eV. Another bright effect of the new physics is 
associated with the isotopic creation of mass by massless fermions (leptons) in graphene. 

Keywords: strong interaction, quarks, gluons, excitons, phonons, QED, QCD 
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introduction. 


The traditional direction of the search for new physics and the origin of mass beyond the 
bound of the Standard Model (SM) in high energy physics is associated with the 
development of more powerful accelerator [1]. A low - energy approach, for example, in the 
problem of the origin of mass, was followed by author of the work [2]. The present brief 
review supplements work [1] with low - energy physics and at the same time is not only 
logical continuation of paper [2] but also its experimental proof [3]. In this connection, we 
must emphasize that outstanding works on hadron mechanics [4] deserve separate and 
fundamental study. 

The history of physics shows that whenever experimental techniques advance to an 
extent that matter, as then known, newer and more powerful studies subsequently show 
that «elementary particles» have a structure themselves (see, e.g. {3 - 6]). Indeed this 
structure may be quite complex, so that the elegant ideas of elementarity must be 
abandoned. The modern concept of the atom emerged at the beginning of the 20th century, 
in particular as result of Rutherford’s experiments [5, 6]. An atom is composed of a dense 
nucleus surrounded by an electron cloud (Fig. 1). 


«, Neutron 
and 


Proton 


Sie 10% 


Atom 


Sires 107? i 


Fig. 1. Structure within the atom. (after http://.Ibl.gov/abs/wallchart/). 


After the discovery of the neutron in 1932 by Chadwick, there was no longer doubt that 
the building blocks of nuclei are protons and neutrons [7, 8] (collectively called nucleons). 
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The electron, neutron and proton were latter joined by fourth particle, the neutrino, which 
was postulated in 1930 by Pauli [6 in order to reconcile the description of B- decay with the 
fundamental laws of conservation of energy , momentum and angular momentum. Thus, by 
the mid - thirties of the 20th century , these four particles could describe all the then known 
phenomena of atomic and nuclear physics. Today, these particles are still considered to be 
the main constituents of matter. in the fifties and sixties was showed that protons and 
neutrons are merely representative of a large family of particles called hadrons [9, 10}. More 
than 100 hadrons have thus far been detected. these hadrons, like atoms, can be classified 
in groups with similar properties. it was therefore assumed that they cannot be understood 
as fundamental constituents of matter. In the late sixties, the quark model (Fig. 1) 
established order in the hadronic zoo [11, 12]. 

Together with our changing conception of elementary particles, our understanding of the 
basic forces of nature and so of the fundamental interaction between elementary particles 
has evolved. As is well - known subatomic physics combines nuclear and particle physics. 
The two fields have many concepts and features in common. The exploration of subatomic 
physics started in 1896 with Becquerel’s discovery of radioactivity; since then it has been a 
constant source of surprises, unexpected phenomena, and fresh insights into the law of 
Nature. In years that followed, the phenomenon was extensively investigated, notably by the 
husband and wife team of Pierre and Marie Curie and by Ernest Rutherford and his 
collaborators, and it was established that there were three distinct types of radiation 
involved. These were named (by Rutherford) @ -, 6 -, and y - rays. We know now that a 
-rays are bound states of two protons and two neutrons (the nuclei of helium atoms), £ - 
rays are electrons and y - rays are high energy photons. The energies of the a -, 6 -, and y - 
particles are of the order of 0.1 MeV up to 10 MeV, whereas energies of the orbital 
electrons are of the order of electron volts. Also, a - particles were found to be barely able to 
penetrate a piece of paper, whereas # - particles could penetrate a few millimeters of 
aluminium, and the y - rays could penetrate several centimeters of lead. 

At about the same time as Becquerel’s discovery, J.J. Tomson in 1897 was first to 
definitely establish the nature of free electrons. the view of the atom at that time was that it 
consist of two components, with positive and negative electric charges, the latter now being 
the electrons. In his experiments on a scattering by gold, Rutherford showed that a - 
particles encountered a very small positively charged central nucleus. To explain the results 
of these experiments Rutherford formulated a planetary model, where the atom was likened 
to a planetary system, with the electrons (the “planets") occupying discrete orbits about a 
central positively charged nucleus ("Sun"). 

Because photons of a definite energy would be emitted when electrons moved from one 
orbit to another, this model could explain the discrete nature of observed electromagnetic 
spectra when excited atoms decayed. In the simplest case of hydrogen, the nucleus is a 
single proton with electric charge +e, where - e is the magnitude of charge on the orbited a 
single electron. Heavier atoms were considered to have nuclei consisting of several protons. 
Examples are carbon and nitrogen, with masses of 12.0 and 14.0. However, it could not 
explain why chlorine has a mass of 35.5 in these units. At about the same time, the concept 
of isotopism was conceived by Soddy. Isotopes are atoms whose nuclei have different 
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masses, but the same charge. The explanation of isotopes was done by Chadwick. He had 
discovered the neutron and in so doing had produced almost final ingredient for 
understanding nuclei. lwanenko and Heisenberg formulated a modern model, according to 
which all nuclei consist of a positively charged proton and a neutral neutron [7, 8]. 

By the early 1930s, the 19th century view of atoms as indivisible elementary particles 
had been replaced and a larger group of physically smaller entities now enjoyed this status; 
electrons, protons and neutrons. To these we should add two electrically neutral particles: 
the photon (7) and the neutrino (v). The photon was postulated by Planck in 1900 to explain 
black ~ body radiation, where the classical description of electromagnetic radiation led to 
results incompatible with experiments. the neutrino was postulated by Fermi in 1930 to 
explain the apparent non - conservation of energy observed in the decay products of some 
unstable nuclei where f - rays are emitted, the so - called f - decay [3]. 

Particle physics is the study of the fundamental constituents of matter and their 
interactions. However, which particles are regarded as fundamental has changed _ with time 
as physicists’ knowledge has improved. The best theory of elementary particles we have at 
present is called, rather prosaically [14], the Standard Model (SM). The SM_ attempts to 
explain all the phenomena of particle physics in terms of the properties and interactions of a 
small number of particles of three distinct types: two spin - 1/2 families of fermions called 
leptons and quarks, and one family of spin -1 bosons - called gauge bosons - ‘which act as 
force carriers’ in the theory. in addition at least one spin-O particle, called Higgs boson, is 
postulated to explain the origin of mass within the theory, since without it all the particles in 
the SM are predicted to have zero mass. All the particles of the SM are assumed to be - 
elementary; they are treated as point particles, without internal structure or excited states. 

Modern physics distinguishes three fundamental properties of atomic nuclei: mass , spin 
(and related magnetic moment), volume (surrounding field strength), which are the source 
of isotope effect. The stable elementary particles (electrons, protons and neutrons) have 
intrinsic properties. Some of these properties such as mass and electrical charge are the 
same for macroscopic objects, Some are purely quantum mechanical and have no 
macroscopic analog. Spin is an intrinsic angular momentum associated with elementary 
particles. The spin angular moment of an electron, measure along any particular direction, 
can only take on the values f/2 or -f/2. The nuclear magnetic moment associates with 
nuclear spin and produces the magnetic interaction with its environment. 

Among the prime premises of quantum physics is the principle of indistInguishabllity of 
elementary particles. All electrons are by definition ‘the same’ and hence, after two (or 
more) electrons interact, it is possible to say ’who is who’. The same applies to other 
elementary particles and even identical atoms (e.g., two Ge atoms are indistinguishable). 
Not so for isotopes. Due to their mass difference, various isotopes of the same chemical 
element (e.g., Ge and ”Ge) are classical distinguishable particles. Within the framework 
of statistical thermodynamics, this classical distingulshability of isotopes renders isotopic 
mixtures a prime illustrative tool for a discussion of the Gibbs paradox. The latter refers to 
the (alleged) discontinuity of entropy upon the mixture of two slightly different species. 

Our present knowledge of physical phenomena suggests that there four types of forces 
between physical bodies (see, e.g. [10, 13}): 
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1) gravitational; 2) electrornagnetic; 3) strong; 4) weak (Table 1). 

Both the gravitational and the electromagnetic forces vary in strength as the inverse 
square of the distance and so able to influence the state of an object even at very large 
distances. Gravitational is important for the existence of stars, galaxes, and planetary 
systems as well as for our daily life, it is of no significance in subatomic physics, being far 
too weak to noticeably influence the interaction between elementary particles [14]. 
Electromagnetism is the force that acts between electrically charged particles (atoms, 
molecules, condensed matter). When nuclear physics developed, two new short - ranged 
forces joined the ranks. These are the nuclear force, which acts between nucleons (proton, 
neutron, etc.) and the weak force, which manifests itself in nuclear f- decay (see, é.g. {10]). 
The nuclear force is a result of the strong force binding quarks to form protons and 
neutrons. Due to experimental results of this paper connected to the manifestation of the 
strong interaction, we should briefly analyze the structure of subatomic particles and the 
strong interaction. We should underlined, that subatomic physics lacks so coherent 
theoretical formulation that would permit us to analyze and interpret all phenomena in a 
fundamental way: atomic physics has such a formulation in quantum electrodynamics, 
which permits calculations of some observable quantities to more than six figures. 
Subatomic physics deals with all entities smaller than the atom. 


Table 1. The four fundamental forces. 


Interaction FQ Mass Range(m) RS Spin TC-S(m?) TTS(s) 
Str Gluon 0 10-5 1 1 10738 10-3 
We W,Z 81,93 GeV/c2 1078 105 1,1 10% 10-8 
Elec Photon 0 00 1/137 1 10°73 10-7 
Gra Graviton 0 re) 10° 2 . ‘ 


Here - FQ field quant, RS relative strength, TC - S Typical cross - section, TTS - Typical 
time scale. 


Particles are identical if they cannot be distinguished one from another by any intrinsic 
property, such as mass, charge, or spin [9]. There are does not exist, in act and in principle, 
any experimental procedure which can identify any one of the particles. In classical 
mechanics, even though all particles in the system may have the same intrinsic properties, 
each may be identified, at least in principle, by its precise trajectory as governed by 
Newton’s laws of motion. This identification is not possible in quantum theory because each 
particle does not possess a trajectory; instead, the wave function gives the probability 
density for finding the particle at each point in space. When a particle is found to be in some 
small region, there is no way of determining either theoretically or experimentally which 
particle it is. Thus all electrons are identical and therefore indistinguishable, as are all 
protons, all neutrons, all hydrogen atoms with'H nuclei, all hydrogen atoms with 7H nuclei, 
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all helium atoms with *He nuclei, all helium atoms with *He nuclei etc. 

Symmetry and its mathematical framework - group theory - play an increasingly 
important role in physics. Both classical and quantum system usually display great 
complexity, but the analysis of their symmetry properties often gives rise to simplifications 
and new Insights which can lead to a deeper understanding. In addition, symmetries 
themselves can point the way toward formulation of a correct physical theory by providing 
constraints and guidelines in an otherwise intractable situation. Taking into account the 
symmetry, the proton and neutron are seen as different state of a single basic nucleon. 
These states bare differentiated by an internal property that have two values , +1/2 and - 
1/2, in analogy with the spin of a particle as the electron. This new property is called isotopic 
spin, or isospin for short, and the nuclear binding force is said to exhibit isospin symmetry 
[12, 13]. Thus, isospin symmetry is an important symmetry in particle physics, although it 
occurs only in the action of the nuclear binding force - or, in modern terminology, the strong 
force. Nuclear supersymmetry is a composite particle phenomenon, linking the properties of 
bosonic and fermionic systems, framed in the context of the Interacting Boson Model of 
nuclear structure [15, 17]. In modern physics, particles are classified into two groups based 
on their properties. There are known as fermions and bosons. In the Standard Model (SM) 
[12, 14] fermions are the fundamental particles of matter. Bosons on the other hand, are 
considered to be the force carriers. Nuclei having an odcl number of nucleons are composite 
fermions whereas nuclei having an even number of nucleons are composite bosons. 

Properties of fermions and bosons are very different especially at temperatures close to 
the absolute zero. Fermions are half - integer spin particles and described by the Fermi - 
Dirac statistics. Fermions obey the Pauli Exclusion Principle. So, two identical fermions do 
not occupy the same quantum state simultaneously. Basically, fermions can be classified 
into two groups: elementary and composite fermions. Elementary fermions are leptons 
(electron, electron neutrino, muon, muon neutrino, tau, and tau neutrino and quarks (up, 
down, top, bottom, strange and charm) (see, Table 2). Hadrons (neutrons, protons) 
containing an odd number of quarks, and nuclei made of an odd number of nucleons (for 
example $C nuclei contain six protons and seven neutrons) are considered to be composite 
fermions. Elementary fermions are the fundamental building blocks of matter and antimatter 
[10, 16]. 


Table 2. Quarks and leptons 


Famil 
ay Electric charge (e) 
1 2 8 
Cr -1 
* Leptons 
Ve Vy Vt 0 


Quarks uc t 273 
ds b -1/3 


- 365 - 


Bosons are identical particles having zero or integer spins. As fermions, bosons can be 
categorized into two groups: elementary and composite bosons. Unlike fermions, bosons do 
not obey the Pauli Exclusion Principle. In other words, any number bosons can occupy the 
same quantum state. Behaviors of bosons are described by the Bose - Einstein statistics. 
The SM [9] only consists of five elementary bosons (Fig. 2). They are namely the Higgs 
boson, gluon, Z and W# bosons [17]. The Higgs boson have zero electric charge and zero 
spin is the only scalar boson [16]. 

The discovery of the neutron by Chadwick in 1932 may be viewed as the birth of the 
strong interaction: it indicated that the nuclei consists of protons and neutrons and hence 
the presence of a force that holds them together, strong enough to counteract the 
electromagnetic repulsion. In 1935, Yukawa [18] pointed out that the nuclear force could be 
generated by the exchange of a hypothetical spinless particle, provided its mass 
intermediate between the masses of proton and electron - a meson. Yukawa predicted the 
pion [9, 10}. The strong forces does not act on leptons (electrons, positrons, muons and 
neutrinos), but only on protons and neutrons (more generally, on baryons and mesons - this 
is the reason for the collective name hadrons). It holds protons and neutrons together to 
form nuclei, and is insignificant at distances greater than 10 ~5 m [10]. Its macroscopic 
manifestations are restricted up to now to radioactivity and the release of nuclear energy. 

The Table 1 given for the strength and range of the forces come from a comparison of 
the effects they produce on two protons. In some respect these resemble an ordinary 
Newtonian force between the protons, varying with the distance between them as if the 
force was derived from a potential function: 

Vin= SF (1) 


for some n. This Is an inverse ~ power force which is diminished by an exponential factor 
at distances larger than a certain distance R, the range of the force. The strength of the 
force is measured by the constant k. The unit of strength is hc/2x% where h is Planck’s 
constant and c the speed of light. Since the protons and neutrons which make up the 
nucleus are themselves considered to be made up of quarks are considered to be held 
together by the color force [11, 12], the strong force between nucleons may be considered 
to be a residual color force (see, also ({10, 14). In the SM, therefore the base exchange is 
the gluon which mediates the forces between quarks. 
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Leptons Quarks 
e, Bt uc, t 
vy, Ys % 


Fig. 2. Summary of interactions between particles described by the Standard 
Model. 


The modern quantummechanical view of the three fundamental forces (all except 
gravity) is that particles of matter (fermions neutrons, protons, electrons) do not directly 
interact with each other, but rather carry a charge, and exchange virtual particles (gauge 
bosons photons, gluons, gravitons) which are the interaction carriers or force mediators. As 
can be see from Table 1, photons are the mediators of the interaction of electric charges 
(protons, electrons, positrons); and gluons are the mediators of the interaction of color 
charges (quarks). In our days, the accepted view is that all matter is made of quarks and 
leptons (see Table 2). As can be see, of the three pairs of quarks and leptons, one pair of 
each ~- the quark u and d and the leptons e and e (velectrons neuino) - are necessary to 
make up the every day world, and a worid which contained only these would seem to be 
quite possible. 

The facts, summarized in the modern nuclear and subatomic physics (see, e.g. [9, 10]) 
allow to draw several conclusions in regard to nuclear forces, most notably that the binding 
energy of a nucleus is proportional to the number of nucleons and that the density of 
nuclear matter is approximately constant. This lead to conclude that nuclear forces have a 
“saturation property". It seems from the last conclusion it is enough to change the number 
of neutrons in nucleus to change strength of nuclear force. But the last one constitutes the 
main ideas of the isotope effect [13]. 


2. Energy bands of electron in a solid. 
2.1. Introduction. 


Any solid consists, i.e. is so collection of nuclei and electrons. In crystalline solids, the 
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atom’s nuclei are located at the nodes of the crystal lattice which has a spatial periodicity. 
The stationary state of all particles is described by the Schddinger equation (see, for 
example [19]) 


Aw=en, (2) 


where A is the Hamiltonian of the entire set of particles, i.e. Hamiltonian of a solid, 
is own wave function, E is the energy of solid body. Let T, T, be the radius-vectors of 


electrons and Ri, R,, .... be the radius - vectors of nuclei. Denote by M, the mass of the 
nucleus of an k - atom and m the mass of the electron. 
The Hamiltonian of the system of particles is described by the following expression: 


fi=R+u (3) 


where K is the operator of the kinetic energy of this system, U is its potential energy. 
The kinetic energy operator for the solid state under consideration will be 


a 
K= ( A; + Xu HA). (4) 


Here A; = & + = + £ js the Laplace operator for i - particle. First term in (4) is the 


operator of ths ‘kinetic send of electrons and the second one of nuclei. The potential 
energy of a set of particles that make up solid is the sum of the energies of the pairwise 
eat of elections with electrons, nuclei with nuclei and electron with nuclei. 


a _Zdie? a ee 
pepo a ra rae 2 Ee 2 cele] ie A. | 2 ud ances |?; - Ae : (5) 


The first two terms in (5) express the energy of the Coulomb repulsion of electrons and 
nuclei, respectively, and the third term is the energy of attracted of electrons to nuclei. 
Substituting (4) and (5) into Schrédinger equation of our system, we get: 


Z,ZLe? = 1 Z,e? 
[4 Fai + La 2m, Ae) +7 2 er FT ate tr a 2 ps | areso [Ae - cal 2 ren ] 
=EY, (6) 
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The wave function in equation (6) depends on the coordinates of all particles, i.e. 


Pe Ti Tas Vases. Pas Pacts, Hav cay Pe (7) 


if this wave function is constrained by its physical meaning (finiteness, unambiguity, 
continuity), then the Schrédinger equation (6) will have a solution not at any energy values 
E, but only at some. These values E, which is solution to equaton (6) determine the 
energies (energy spectrum) of the solid. However, due to the vast number of independent 
variables, equation (6) cannot currently be solved in general form. To find an approximate 
solution, they resort to a number of simplifying assumptions. 

Given the large difference in the masses of nuclei and electrons (M, >> m) and therefore 
the nature of movement these particles is significantly different. This approximation is called 
the adiabatic or Born - Oppengeimer {20] approximation. indeed, if the nuclei rest, then the 
kinetic energy of the nuclei turns to zero. The potential energy of interaction of nuclei 
becomes some constant, i.e. 


1 228? aa 
py py “ao -F] const. (8) 


By choosing the starting point, we can turn it to Zero. Taking this into account, the 
Schrédinger equation accepts the form: 


& A; Sk 28? __. = 
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(9) 


Here the energy of electrons E, and their wave function ‘Y, depend on resting nuclei Rox 
only parametrically. The following approximation is called valence approximation 

Despite significant simplifications, the Schrédinger equation (9) cannot be solved. The 
following approximation is called valence approximation. It is believed here that all electrons 
of the inner shells of an atom form a resting atomic residue together with nucleus, and, 
therefore, the Schrédinger equation is written only for valence electrons. As part of the 
adiabatic approximation and valence approximation, the wave function of the system 
remains dependent on the coordinates of a valence electrons. Since the latter interact with 
each other, the variables in the Schrédinger equation (9) are not separated . Therefore, a 
single - electron approximation is used to solve the problem, for which the Hartree - Fock 
method is used. The main ideas of this method is to replace the potential energy of 
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which is the interaction energy of the i - electron with some effective field in which each 
electron moves independently. This effective field characterizes the effect of all other 
electrons on the i - electron. Indirectly, it also depends on the movement of the i - electron 
affects the movement of all other electrons. Assuming that such a field is found, rewrite the 
Schrédinger equation (9) as 


{x[- EA; + OP) + UP) J}. =E.Y.. (10) 


i 


Here the U,(7;) denotes the potential energy of the i - electron in the field of all nuclei. 
Under the sign of the sum (10) is the Hamiltonian i - electron 


A; =- 4A; + Ci?) + U,T). (11) 
So the Schrédinger equation can be rewriten like this: 


AY, =(SA)e. =E,W, (12) 


Since now Hamiltonian does not contain the energy of interaction of electrons and is the 
Hamiltonian sum of individual electrons, the solutions to equation (12) is the product of 
single - electron functions 


WAP, Pa, o.) =wi(Pi wo(Pr).... =Myi(P:). (13) 


And each function wi(Ti) satisfies the one-electron Schrédinger equation 
Av: =Ewi, (14) 


in which the interaction of the i - electron with the rest is described by the potential 
Gi) 


Thus, the introduction of an effective field allows to reduce the multi - electron equation 
to a single - electron system. At the same time, the energy of systems 


E, = > Er. (15) 
i 
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Although the wave function (13) is a solution to the Schrodinger equation for a crystal, it 
does not satisfy The Pauli principle [21]. According to Pauli’s principle, one quantum state 
cannot contain more than two electrons with different spin orientations. The complete wave 
function of the system satisfying this condition must be antisymmetric, i.e. when rearranging 
coordinates and projecting the spin, it must change the sign. Our function (13) does not 
satisfy this condition. As a rule, an antisymmetric function is written in the form of a Slater 
determinant (for more details see, for example (21]). The antisymmetric properties of such a 
function derive from the properties of the Slater determinant. Since atom are located in 
space strictly periodically in the crystal, the full potential V(?’) = ured ) + U(T) must have a 
three - dimensional periodicity. Although the exact species V(7) is unknown, it is important 
to know that this potential is a periodic function whose period coincides with the period of 
the crystal lattice of the object under study. This periodicity leads to the formation of a band 
energy structure (Fig. 3). An electron from the valence band (see Fig. 3) is excited into the 
conduction band. 


aT 
EEL 


QO wavector n/a space coordinate x., density of states D 


Fig. 3. Various possibilities to present the band - structure of homogeneous, 
undoped insulator (semiconductor). 1 - the dispersion relation, i.e; the energy E as a 
function of the wave vector K, 2 - the energy regions of allowed and forbidden states 
as function of a space coordinate x and , 3 - the density of states (all curves are 
schematic ones). 


The dependence of the energy of an electron on its quasi - pulse (the curvature of 
isoenergy surfaces in rae space) is connected by another fundamentally important 
characteristic of the electron, which is called its effective mass. According its definition [19] 
we have: 
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The effective mass reflects the influence of the lattice potential on the motion of the 
electron in the crystal under the action of an external force. In thé case where the electron 
is in the vicinity of the minimum energy, i.e. the vicinity of the bottom of the conduction band 
a ) 0 and m* = constant ) Q, i.e. the electrons behave as negatively charged particles 


with a positive mass. In the vicinity of the maximum energy, i.e. in the vicinity of the vertex 
of the valence band oe. ( 0 and m* = constant { 0 and, therefore in the vicinity of the 
dk 


vertex of the valence band , such a charge carrier behaves as a particle with a positive 
charge and a positive effective mass and is called a hole. As will be shown below by way of 
specific examples, the effective electron mass in different crystal varies from 1 to hundreths 
of the free electron mass. 


2.2. Excitons. 


The properties of the electron and the hole are both described by the band structure 
within the one - electron approximation. in this section we shall go beyond this 
approximation and consider the effects of electron - electron interaction on the absorption 
spectra. To simplify the calculation we shall make the following assumptions. We shall 
include only the Coulombic part of the electron ~ electron interaction neglecting both 
exchange and correlation terms. Furthermore, the interaction between the excited electron 
in conduction band (see, Fig. 3) and those left behind in the now almost filled valence band 
will be replaced by an electron - hole interaction. Attraction between the electron and the 
hole causes their to be correlated and the resultant electron - hole pair is known as an 
exciton. 

It has been more than eight decades since the introduction of quasiparticle exciton by 
Frenkel [22] and the extreme fertility of this idea has been demonstrated most powerfully. 
According to Frenkel the exciton is an electron excitation of one of the atoms (ions) of the 
crystal lattice, because of the translation symmetry, moves through the crystal in an 
electrically neutral formation. Since Frenkel the concept of an exciton has been developed 
in the papers of Peierls [23] and Slater and Schokley [24]. Problems concerning light 
absorption by solid state have been considered somewhat differently Wannier [24] and Mott 
[25]. According to the Wannier - Mott results the exciton is the state of an electron and hole 
bonded by the Coulomb force. The electron and hole in exciton state are spatially separated 
and their charges are screened. In the Frenkel papers the excitations localized on the lattice 
site were described thus, after the Wannier - Mott papers, the excitons became divided into 
the excitons on the Frenkel (small radius) excitons (for details see [55]) and the Wannier - 
Mott (large radius) excitons ( [65]). However, a description of the basic difference between 
these two models is absent (27 - 29]. The experimental discovery (see e.g. Gross [30]) of 
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the Wannier - Mott exciton (see Fig.4) on the hydrogen-like absorption spectrum in the 
semiconducting crystals was the basis of a new subject - exciton physics (see also Agekyan 
[31]; Permogorov [32]). The influence of external perturbation (electrical and magnetic 
fields, uniaxial and hydrostatic deformation) on the optical spectra of the Wannier - Mott 
excitons (see e.g., Gross [30]) and their energetic characteristics has been demonstrated 
repeatedly. These investigations permitted high-accuracy measurements not only the 
exciton binding energy but also of their translational mass, values of effective masses of the 
electron and hole, their g - factors etc. Moreover, the detailed account of the photon - 
exciton interaction has led to the concept of polaritons [33]. From the time of the 
experimental discovery of the Wannier - Mott exciton the problem concerning the interaction 
of excitons and the crystal lattice has persisted for more than four decades (Haken [34]; 
Haug and Koch [35)). 


Fig. 4. Discrete and continuos (hatched area) Wannier - Mott exciton energy 
spectrum taking into account its kinetic energy. The broken line connects to the 
dispersion o light in the medium. R - is the exciton binding energy, n = 1, 2, 3, ..... 


Below we briefly consider quantitatively the effect of Coulomb attraction on the motion of 
electrons and holes (see, also [36]) in the vicinity of an Mo critical point of a direct bandgap 
semiconductor in three dimensions. Further we shall assume the conduction band to be 
spherical with energy 
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E.(m) =e, + M2, (17) 


where E, is the band gap, and the corresponding hole energy to be given by 


E,(k)-E toe. (18) 


Let the Bloch function and the hole be represented by ¥> (7¢) and Wp 


(7h n)s respectively. As usually we assume that the Coulomb interaction between slacken 
and hole is weak due to screening by the valence electrons so that the effective mass 
approximation is valid. We can write the exciton wavefunction ‘¥ as linear combination of 
the electron and hole wave functions: 


v= 0 o(k i 1) (73)¥z (7). (19) 


Ke, kn 


Similar of the donor electron (see, e.g. [37, 38]), the electron and hole in an exciton are 
localized relative to their center of mass, so it is more convenient to express their 
wavefunctions in terms of Wannier functions rather than Bloch functions. in terms of the 
Wannier functions a;- (re) and ane (ff) for electron and hole, respectively, the exciton 


wavefunction can be written as 


(Fe fh) = ND) (Re, Rn) ag: (Fe)ag (TR), (20) 


Re. Rh 


where (Re, Rn) is the exciton envelope wavefunction. The wave equation for 
(Re, Rh) is 


[- Bia ) Vee * (ais Ve —$- |o(Fe, Fn) =E0(Re, An) , (21) 
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where «9 is the zero ~- frequency dielectric constant of material (insulators or 
semiconductors). Equation (21) can be solved in the same way as in the case of the 
hydrogen atom [39]. 

One expresses Re and Ry, in terms of two new coordinates: a center of mass coordinate 
R and a relative coordinate?’ defined by 


— = 
_ ReMme +Ramn > RL] 
R- ae and ? = Re - Rn. (22) 


The equation of motion for the center of mass is now decoupled from that for the relative 


motion because the Coulomb interaction term does not involve R (see, also [36, 38]). The 
two resultant equations are 


(-47)v3¥(F) = Env(R), (23) 
(-E-V2 - £)6(?) =End(?), (24) 


where y:, the reduced mass of the exciton, is defined by 
+ = aig + wt: (25) 


The total energy of the exciton E is simply the sum of Er and E;. The solution of (23) 
and (24) can be obtained readely. Equation (23) describes a free particle whose 
eigenfunction and energy are given by 


Wk (R) =N*exp(iK A) and Er= 44. (26) 


Er represents the kinetic energy of the center of mass motion [86]. 

Equation (25) is similar to the equation describing the motion of the donor electron [38]. 
As in the hydrogen atom, its wavefunctions and energies can be described by three 
quantum numbers: a principal quantum number n, the angular momentum quantum number 
! and the magnetic quantum number m. The wavefunction ¢ can be expressed _in polar 
coordinates (r, 8, @) as 
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daim(?) =Rni{?)¥im(8, 9), (27) 


where'Rni(?) and Yim(8, ¢) are, respectively, the associated Laguerre polynomials 
and the spherical harmonic functions. These functions are tabulated in many quantum 
mechanics textbooks and therefore will not reproduce here (see, e.g. [ 40]). For isotropic 
effective masses E,; depends on n only and is given by 


E,(n) =E,(co) - &., (28) 


where E;(o) is the minimum energy of the continuum states (see, Fig. 4), i.e. the 
energy gap Eg, and R* is the Rydberg constant for the exciton defined as [36] 


R* = fry = 13.6 (+4) ev. (29) 


Combining the above result for the relative motion and the center of mass rnotion of the 
exciton, we obtain the following envelope functions and energies for the excitons: 


Orim(A, P) = (1//M)exp(iR «RR (7)¥m(8, ), (30) 


Enim = Eg + oe Be, (31) 


The energy spectrurn of a Wannier - Mott excitons is shown in greater detail in Fig. 4. 
The above model of excitons based on electrons and holes with spherically symmetric 
parabolic dispersion is useful for understanding exciton effects on optical spectra. 
However, it is not accurate enough for quantitative interpretation of experimental spectra in- 
diamond and zinc - blende - type semiconductors. Of the various atternpts to calculate 
excitonic effects based on realistic band structures, we shall mention the paper [41] (see, 
also [42)). 


A) LiH 


LiH crystal is of NaCl structure type. The simple electronic structure of Lit and H~ ions, 
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having 1s? configuration, gives LiH a special place among the binary crystals and in many 
aspects allows it to serve as an ideal model for ionic compounds. The spectrum of one - 
electron states of crystal is determined by the solution of Schrodinger equation for the 
‘extra’ particle (hole or electron) moving in the averaged field created by all the remaining 
electrons and nuclei (see, also [43]): 


[-A + VE) J ¥e@) = En (KYM, (32) 


where the notation is conventional, and V(r) is the periodic potential. The existing 
calculation techniques differ in the method of constructing the electron potential V(7), the 
approximation of the wave function ‘¥;(F), the ways of ensuring self - consistency, the 


reliance on empirical parameters, etc. In particular, two factors are especially important in 
case of ionic crystals: (1) the inclusion of exchange interaction, and (2) the inclusion of 
polarization of the electron and ion subsystems of crystal by the extra particle [44]. 

The one-electron potential of any many-electron system is nonlocal because of the 
exchange interaction between the electrons. It is very difficult to take this interaction into 
account. Because of this, the exact potential in the band theory is often replaced by the 
local potential of the form Veschange(7) oc a{p(*)}'", where p(7) is the charge density function, 
and the constant a is selected in the range from 1 (Slater potential [45]) to 2/3 (Cohn-Sham 
potential) [43] 

Two effects are associated with the local exchange. First, the results of calculation 
depend strongly on the numerical value of a, and second, this approximation always 
underestimates the values of E, and E, (the width of the valence band). It is the low values 
of E, obtained by many authors in the approximation of local exchange that are responsible 
for the wrong conclusion concerning the inapplicability of band theory to tonic dielectrics. 
These problems reflect the fundamental drawback of the one-electron approximation which 
does not take into account the reciprocal effect of the selected electron (hole) on the rest of 
the system. This effect consists in the polarization of the crystal by the particle, and is 
generally made up of two parts: the electron polarization (Inertialess), and the lattice 
polarization (inertial). In the common optical phenomena, related to absorption or scattering 
of photons, the lattice polarization is not important, because the frequency of optical 
transitions is much higher than the average frequencies of phonons. The electron 
polarization is different. The extra particle (electron, hole) is regarded by this theory as the 
slowest particle in the system - in other words, all the remaining electrons adiabatically 
follow it. Hence it follows that the inertlaless polarization definitely must be included in the 
calculation of energy spectrum. An important feature of ionic crystals is that the polarization 
energy E, is of the same order of magnitude as the bandwidth. Such a correction obviously 
cannot be regarded as small. In the extreme case of particle at rest, the polarization energy 
can be calculated by methods of classical electrostatics (the Mott - Littleton method [19}), or 
by the newer and more accurate technique proposed by Fowler [46]. The value of E, for 
AHC found by this method is 2-3 eV for each of the quasi-particles (E, > 0 for electrons, 
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and E, < 9 for holes). This implies that the inclusion of electron polarization will reduce the 
magnitude of E, by 4-5 eV [43]. By assumption, the electron bands are displaced rigidly, 
without changing the dispersion law E(z). 

The simple electron structure of lithium hydride (combined with the negligibly small 
spin-orbital interaction) is very helpful for calculating the band structure: all electron shells 
can easily be taken into account in the construction of the electron potential. The first 
calculations of band structure of lithium hydride were carried out as early as 1936 by Ewing 
and Seitz [47] using the Wigner-Seitz cell method. This method consists essentially in the 
following. The straight lattice is divided into polyhedra in such a way that the latter fill the 
entire space; inside each polyhedron is an atom forming the basis of the lattice 
(Wigner-Seitz cells), The potential inside each cell is assumed to be spherically symmetrical 
and coinciding with the potential of free ion. This approximation works well for ions with 
closed shells. The radial Schrodinger equation in the coordinate function R,(r) is solved 
within each selected cell, the energy being regarded as a parameter. Then the Bloch 
function is constructed in the form of expansion 


WEP) = Doe Dorp Cm (K) Yin (8, 9) RIC? E), (33) 


where 7, 8, @ are the spherical coordinates (with respect to the center of the cell); Yi», are 
spherical functions. The coefficients Cin(k) and the energy E(k) are found from conditions of 
periodicity and continuity on the boundaries of the cell. if 7; and 7% are the coordinates of 
two points on the surface of Wigner-Seitz cell, linked by the translation vector Ri, then the 
boundary conditions are 


V2) = exp(ikRi) ¥ 2), Sa ae 
Vi e72) = exp (-ikR, V2), (35) 


where V,, is the gradient normal to the surface of the cell. We see that the method of 
cells only differs from the problem of free atom in the boundary conditions. Owing to the 
complex shape of the cell, however, the construction of boundary conditions is a very 
complicated task, and this method is rarely used nowadays. 

The method of plane associated waves (PAW) was used for calculating the band 
structure and the equation of state for LiH was used in Perrot [48]. According to this 
method, the crystal potential is assumed to be spherically symmetrical within a sphere of 
radius 7; described around each atom, and constant between the spheres (the so-called 
cellular muffin-tin (MT) potential). inside each sphere, like in the Wigner-Seitz method, the 
solutions of Schrodinger equation have the form of spherical harmonics; outside the 
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spheres they become plane waves. Accordingly, the basis functions have the form 
We) = exp (ik?) 0(F - 7.) + CamYm(8, 9) RCE, 7) OCF; -#), (36) 


where 6(x) = 1 at x > 0, and 6(x) = 0 at x < 0. The coefficients a,,, can be easily found 
from condition of sewing on the boundary of the sphere. This is an important advantage of 
the PAW method over the method of cells. The calculations of Perrot [48] are 
self-consistent, and the local potential is used in the Cohn-Sham form. The correlation 
corrections were neglected. The method of Corringi-Cohn-Rostocker (CCR method), or the 
method of Green’s functions, was used for calculating the band structure of LiH in Zavt et 
all. [49] (only concerned with the valence band) and in Kulikov [50]. Calculation of band 
structure of LiH in Grosso and Paravicini [51] was based on the wave function used in the 
method of orthogonalized plane waves (OPW) of the form 


WG) = exp(ik?) - 0, (exp(ik?) | X,)X@), (87) 


where X. are the atomic functions of state of the skeleton; (exp(ik?) | X,) is the 


integral of overlapping of plane wave with skeleton function (see also [52]). The method of 
linear combination of local basis functions was applied to the calculation of band structure of 
LIH in Kunz and Mickish [53]. This method is based on constructing the local orbitals for the 
occupied atom states, based on certain invariant properties of the Fock operator. The main 
feature of local orbitals is that they are much less extensive than the atom orbitals. 
Importantly, the correlation correction is taken into account in Kunz and Mickish [53]. Owing 
to the high polarizability of hydrogen molecules, the correlation effect in lithium hydride is 
exceptionally strong. Yet another calculation of band structure of LiH was carried out in Zavt 
et al. [49] using the so-called method of extended elementary cell. This approach is based 
on the semlempirical techniques of the theory of molecules, and is much similar to the 
cluster calculations. Let us add that the cluster is selected in such a way that the 
quasi-molecular wave function transforms in accordance with the group symmetry of certain 
wave vectors in the Brillouin zone, This methods only yields the energy values at the points 
of high symmetry. We ought to mention also Hama and Kawakaml [54], where, in 
connection with the study of high pressure effects on the transition NaCl - CsCl! in lithium 
hydride, the band structure and the equation of state of the latter are analyzed in detail. The 
calculated band structures of LiH are compared in Fig. 5. 


Fig. 5. band structure of LiH crysta als calculated: [53] - 1; By - 2; [52] - 3. 


We see that the overall picture given by various methods is generally the same, despite 
the vast spread of the transition energy values (see Table 3). 


Table 3. Calculated energy values of some direct optical transitions in LiH 
reduced to the experimental value of E, = 5.0 eV. 


Transition 1 2 34 5 
K, - K3 6.9 7.5 65 6.4 
W:-W3; 80 7.9 7.3 7.4 
Li-L2 92 96 9.0 9.1 
Wi-W. 126 14.9 12.2 

X,- X's 12.9 13.8 13.6 

Ki ~ Ka 14.7 16.1 15.0 
Li-L’s 19.7 20.9 20.7 
T,-Tis 245 25.3 33.3 


Looking at the structure of the valence band we see that it is very similar to the s - band 
in the method of strong bond (see also [87}). This is surprising, given the strong overlap of 
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the anion s - functions in lithium hydride. The wave functions in this band are almost entirely 
composed of the Is states of hydrogen ion. Different authors place the ceiling of the band 
either at point X or at point W of the Brillouin zone. Although in all cases the energies of the 
states X; and W; differ little (< 0.3 eV), the question of the actual location of the top of the 
valence band may be important for the dynamics of the hole. Different calculations also 
disagree on the width of the valence band. For example, the width of the valence band in 
LiH without correlation is, according to Kunz and Micklsh [53], Ey = 14.5 eV, and the value 
of Ey, is reduced to one half of this when correlation is taken into account. This shows how 
much the polarization of crystal by the hole affects the width of the valence band Ey. 
According to Perrot [29], the width of the valence band in LiH is 5.6 eV. The density of 
electron states in the yalence band of LiH was measured in Betenekova et al. [55] and 
Ichikava et al. [56]. In Betenekova et al. the measurements were carried out with a magnetic 
spectrometer having the resolution of 1.5 eV, whereas the resolution of hemispherical 
analyzer used in Ichikawa et al. [56] was 1.1 eV. From experimental data, the width of the 
valence band is 6 eV according to Betenekova et al., and 6.3 eV according to Ichikawa et al. 
Observe the good agreement with the calculated value of Ey in this theory. Let us add also 
that the measured distribution of the electron density of states in the valence band of LiH 
exhibits asymmetry typical of s-bands (for more details see Betenekova et al. [55] and 
Ichikawa et al. [56]). The lower part of the conduction band is formed wholly by p-states and 
displays an absolute minimum at point X which corresponds to the singlet symmetry state 
X4. The inversion of order of s and p - states in the spirit of LCAO method may be 
understood as the result of the s - nature of valence band. Mixing of s - states of the two 
bands leads to their hybridization and spreading, which changes the sequence of levels. If 
we compare the structure of the conduction band with the p - band of the method of strong 
bond (see, also [58]), we see that the general structure and the sequence of levels are the 
same except for some minor details (the location of Ls level, and the behavior of E(k) in: the 
neighborhood of T1s). In other words, the lower part of the conduction band in lithium 
hydride is very close to the valence p - band of AHC. The direct optical gap In LiH according 
to all calculations is located at X point and corresponds to the allowed transition X;-X4. The 
indirect transition W1-X4 ought to have a similar energy. According to the above 
calculations, the energies of these transitions differ by 0.03-0.3 eV. The different values of 
E, for LiH obtained by different authors are apparently due to the various methods used for 
taking into account the exchange and correlation corrections (see above). As follows from 
Table 3, the transitions at critical points in the low-energy region form two groups at 7 - 9 
and 13 - 15 eV. Measurements of reflection spectra in the 4 - 25 eV range at 5 K (Kink etal. 
[57]) and 4 - 40 eV at 300 K throw new light on the results of calculations (see also review 
by Plekhanov [58]). The singularities occurring at 7.9 and 12.7 eV in reflection spectra are 
associated in the above papers with the interband transitions W, -Wa and X: -Xs 
respectively. ; 

From the standpoint of dynamics of quasi-particles, an important consequence of such 
band structure is the high anisotropy of the tensor of effective mass of electrons and 
(especially) holes. The estimated mass of electron in the neighborhood of X4 is, according 
to Kunz and Mickish [53], (7), = 0.3mMo in the direction X - P, and (m.)y = (m_)z = 0.8 Mo in 
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the direction X - W. Similarly, the mass of hole in the neighborhood of X, is X - Tin the 
direction (Me)x = 0.55 Mo and about the same in the neighborhood of W;. It is assumed that 
the transverse components of m, are greater by several orders of magnitude (Zavt et al. 
[49]). Note also that, according to Baroni et al. [52], the estimated masses of carriers are: 
Met = 0,121; Me = 0,938; My:= 0,150; my= 4,304 me, where the subscripts 1 and t denote, 
respectively, the longitudinal (in the direction r - X) and the transverse (in the direction X - 
W) components. This high anisotropy of masses of electron and hole ought to have resulted 
in the high anisotropy of the reduced (1/j: = f/m. + 1/m;) and the translation (M =m, + my) 
masses of exciton. This, however, is not the case. Moreover, the study of Plekhanov and 
Altukhov [59] reveals that with a good degree of confidence one may assume that in the 
energy range E < 40 meV the exciton band is isotropic and exhibits parabolic dispersion (Me 
= 0.04m and m, = 0.15m). As was shown below, the studies of exciton ~ phonon 
luminescence of free excitons and resonance Raman scattering of light in LiH crystals [43] 
reveal that the kinetic energy of excitons in these crystals is greater than E, by an order of 
magnitude exactly because of the very small masses of electron and hole. The last one 
may indicate that in the metallic phase of hydrogen at the high pressure [43] we can expect 
the Dirac character of the electronic excitations [60]. We should add that isotope 
substitution will be very useful method for renormalization of the band - gap energy in 
graphene - future semiconducting material [61]. 


B) Diamond 


Carbon atom is built from 6 protons, A neutrons and 6 electrons, where A = 6 or 7, yield 
the stable isotopes '2C and C, respectively, and A =8 characterizes the radioactive isotope 
4C. The isotope 'C, with nuclear spin | = 0, is the most common one in nature with 99% of 
all carbon atoms, whereas only ~ 1% are #C with nuclear spin | = 1/2. There only traces of 
'4C (107!? of all carbon atoms) which £ - decays into nitrogen '*N. Although '4C only occurs 
rarely, it is important isotope used for historical dating (see, e.g. [62]). 

Carbon, one of the most basic elements in nature, still gives a lot surprises. It is found in 
many different forms - allotropes - from zero dimensional fullerene, one dimensional carbon 
nanotubes, two dimensional graphene and graphite, to three dimensional diamond (Fig. 6) - 
and the properties of the various carbon allotropes can vary widely [63]. For instance, 
diamond is the hardest material, while graphite is one of the softest: diamond is transparent 
to the visible part of spectrum, while graphite is opaque; diamond is an electrical insulator, 
while graphite is a conductor. Very important is that all these different properties originate 
from the same carbon atoms, simply with different arrangements of the atomic structure. 
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Fig. 6.Structure of some representative carbon allotropes. 


The Diamond is one of the simplest and most important systems, especially in isotope 
effect of condensed matter. It has considerable mechanical strength; also it has high 
thermal conductivity and chemical stability even at several hundreds of degrees Celsius (64, 
65]. Diamond has, perhaps, the simplest and most basic covalent band structure. Fig. 7 
shows the electronic energy - band along several lines of high symmetry from the center (I) 
to the boundary of the first BZ as calculated with modified LCAO method in paper [66]. 


Oba vector k 


- 383 - 


Fig.7. Electronic energy band structure of diamond along several lines of high 
symmetry from the center to the boundary of the first BZ as calculated with a 
modified LCAO method [66]. 


The electronic states are labeled using the notation for the single group of the diamond 
structure. The location of several interband transitions are included by the vertical arrows. 
The fundamental absorption edge of diamond corresponds to indirect transitions from 
highest valence band at the © point to the lowest conduction band in the A direction (X 
point) (i.e., P25: -> Ai [X;,]). The theoretical indirect gap energy E” lays between 5.05 and 
5.6 eV [67], that in reasonable agreement with the experimental data: E/? = 5.470 + 0.005 
eV [64, 65]. LIke in Si, the lowest - lying conduction band at© point of BZ in diamond, M15 is 
p - like, however, as we can see, in Ge the s - like I’. band is the lowest conduction band. 
As mentioned in literature (see, e.g. [68] and references therein), one common 
characteristic of all published spectra in the UV region of diamond is the presence of two 
major peaks at ~7.02 - 7.4 eV and ~ 12.2 - 12.7 eV. The first peak may originate from I'2.: 
-~» Tys transition (E»’) and the second one from X4 — X; and &2 — &s transitions (E2) (the 
more detail see [65)). 

The dependence of the band gap energy on isotopic composition has already been 
observed for insulators and lowest (indirect-direct) gap of different semiconductors ([69, 
70]). It has been shown to result primarily from the effect of the average isotopic mass on 
the electron-phonon interaction (for details see below), with a smaller contribution from the 
change in lattice constant. This simplest approximation, in which crystals of mixed isotopic 
composition are treated as crystals of identical atoms having the average isotopic mass is 
referred to as virtual crystal approximation (VCA). Going beyond the VCA, in 
Isotoplically mixed crystals one would also expect local fluctuations in the band gap energy 
from statistical fluctuations in local isotopic composition within some effectlve volume, such 
as that of an exciton. 


C) Graphene. 


Graphene is a single - atom - thick two - dimensional planar layer of carbon atoms in 
hexagonal honeycomb array composed of two superposed triangular sub - lattices. 
Graphene, a monolayer of sp? - bonded carbon atoms or one monolayer of graphite. It can 
be wrapped up into OD fullerenes, rolled into 1D nanotubes or stacked 3D graphite. The 
band structure of graphene involves two nodal zero - gap points (K - points, see Fig. 8) in 
the first Brillouin zone at which the conduction and valence band touch. In the vicinity of 
these points, the low - energy electron/hole energy dispersion relation is proportional to 
momentum, rather than its square. This is analogous to the energy dispersion relation of 
massless relativistic electrons, so the electrons/holes of graphene are described as “Dirac 
fermions" having no mass. The relativistic - like description of electron waves on 
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honeycomb lattices has been known theoretically for many years, never failing to attract 
attention, and the experimental discovery of graphene now provides a way to probe 
quantum electrodynamics (QED) phenomena by measuring graphene’s electronic 
properties (for details see [15] and references therein). 

Graphene can be considered as the building block of many carbon allotropes. As was 
shown [15] (Fig. 3°, chapter 5 in [15]) it is a two - dimensional crystal with hexagonal 
structure consisting of bipartite lattice of two triangular sublattices. We should remind that 
the fourth valence electron Is in the 2p, orbital which is orthogonal to the graphene plane. A 
weak z bond is formed by the overlap of half - filled 2p, orbitals with other 2p, orbitals. The 
transport properties of graphene are determined by these delocalized a electrons [71]. 
Further we will present a very brief description of the band ¢ructure. As mentioned above, 
graphene is a honeycomb lattice of carbon atom (Fig. 6). it is a bipartite lattice with two 
sublattices A and B that are triangular Bravais lattice with two atoms per unit cell (A and B). 
Each A - or B - type atom is surrounded by three atoms of the opposite type. In a simple 
neighbor model graphene is a semimetal with zero - overlap between valence and 
conduction bands [15]. Considering only the xy plane, the unt vectors in real space, @: and 
@2, and the reciprocal lattice vectors Bi and b are shown in (Fig. 3, chapter 5 in [15]. The 
real space lattice vectors are written as 


a= 2(J%, 1), a, =4(J/3, -1), (38) 


where a= | a; | = | @2 | = 0.246 nmis the lattice constant. The corresponding unit in 
the reciprocal lattice are 


Di= (+, 1), B= (4, 1), (39) 


with a reciprocal space lattice constant 4z/(aJ3 ). The first Brillouin zone is a hexagon 
(Fig. 3* in chapter 5 in [15]), where the corners are called the K points (see below). 

The energy dispersion of z electrons in graphene was first derived in 1947 by Wallace 
[72] within the tight - binding approximation. In this case, the wave function of graphene is a 


linear combination of Bloch function for sublattice A 


O4= Deol? - Ra), (40) 
Ra 
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and equilibrium function ©, for the B sublattice. Here N is the number of unit cells, RA, 
are the position of the atom A and 9(? - Ra) is the 2p, orbital of the atom at A. The sum 
runs over all unit cells, i.e. all possible lattice vectors. in the nearest neighbor approximation 
(every A site has three nearest B sites, and vice versa), the energy eigenvalues can be 
obtained in a closed form [72, 73] 


e(kx, Ky) = + ¥o [1 + 4cos-2 4 cos“ + 4cos?-4112, (41) 


where yo is the transfer integral between the nearest neighbors. The energy dispersion 
of two - dimensional graphene according to this formula is plotted in Fig. 8 as a function of 
the wave vector K. 


Fig. 8 Energy dispersion of graphene obtained within the tight - binding 
approximation. (a) Energy dispersion relation for graphene, drawn in the energy 
region of the Brillouin zone. Since in this approximation to ignore the coupling. 
between the graphene sheets, the band depend only on k, and k,. The z band is 
completely filled and meets the totally empty <* band at the K points. Near these 
points both bands have linear dispersion as described in the text. (b) The dispersion 
along the high - symmetry points MK. 


The upper half of the curves is called the z* or the antlbonding band while the lower one 
is z or the bonding band. The two bands degenerate at the two K points given by the 
reciprocal space vectors K = (2m/a)(1/3, 1/3 ) and K’= (2n/a)(- 1/3, 1//3 ) points where 
the dispersion vanishes (Fig. 9). 


Fig. 9. Dirac cones at the K points. 


Three o —bonds join a C atom to its three neighbors. They are quite strong, leading to 
optical - phonon frequencies much higher than observed in diamond (see below). Such 
orbitals are even with respect to the planar symmetry. The s bonds are strongly covalent 
bonds determining the energetic stability and the elastic properties of graphene. The 
remaining p, orbital, pointing out of the graphene sheet is odd with respect to the planar 
symmetry and decoupled from the s states. From the lateral interaction with neighboring p, 
orbitals (called the ppz interaction), localized x (bonding) and x* (antibonding) orbitals are 
formed [73]. Graphite consists of a stack of many graphene layers. The unit cell in graphite 
can be primarily defined using two graphene layers translated from each other by a C-C 
distance (a... =1.42 A). The three-dimensional structure of graphite is maintained by the 
weak interlayer van der Waals interaction between x bonds of adjacent layers, which 
generate a weak but finite out-of-plane delocalization [73]. The bonding and antibonding c 
bands are actually strongly separated in energy > 12 eV at , and therefore their contribution 
to electronic properties is commonly disregarded, while the bonding and antibonding x 
states lie in the vicinity of the Fermi level (see, e.g. Fig. 2 in [74]) . The two remaining x 
bands completely describe the low energy electronic excitations in both graphene and 
graphite (see [73] and references therein), The bonding x and antibonding x* orbitals 
produce valence and conduction bands which cross at the charge neutrality point (Fermi 
level of undoped graphene) at vertices of the hexagonal Brillouin zone. Carbon atoms in a 
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graphene plane are located at the vertices of a hexagonal lattice. 

Basically, graphene has redefined the limits of what a material can do: it boasts record 
thermal conductivity and the highest current density at room temperature ever measured (a 
million times that of copper!); it is the strongest materlal known (a hundred times stronger 
than steel!) yet is highly mechanically flexible; it is the least permeable material known (not 
even helium atoms can pass through it!); the best transparent conductive film; the thinnest 
material known; and the list goes on [75]. in the vicinity of K - points of the Brillouin zone, 
the low - energy electron/hole dispersion relation is proportional to momentum, rather than 
its square. This is analogous to the energy dispersion relation of massless relativistic 
electrons, so the electrons/holes of graphene are described as Dirac fermions having no 
mass (see also Fig. 10). in a simple neighbor model graphene is a semimetal with zero - 
overlap between valence and conduction bands. in order to make graphene a real 
technology, a special issue must be solved: creating an energy gap at K - points in the 
Brillouin zone (63, 75, 76]. 
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Fig. 10. Cyclotron mass m, of electrons and holes as a function of their 
concentrations Symbols are experimental data, solid curves the best fit to theory 
(after [77]). 
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Further we describe how the tight - binding model discussed above yields the famous 
massless Dirac fermions which have made the study of monolayer graphene so enticing. 
As noted above, the Fermi energy corresponds to « = 0 at the K points. The low - energy 
properties, corresponding to the electronic states near the Fermi energy, can be described 
by expanding the energy dispersion around the K points. Writing the graphene wave vector 
q =K + K, such that | K |a« 1, one can write a Taylor expansion of :(q) near K which 
yields 


St (K) = +8 yoak = tVehk, (42) 


where vr = /3yoa/2h x 10° ms~ is the Fermi velocity. This velocity is independent of 
the carrier density {77] therefore the energy dispersion corresponds to two cones meeting 
at the K point (see Fig. 8%) with linear dependence on the wave vector. The corresponding 
density of states is 


ple) «] e |, (42) 


i.e. the density of states is linear and vanishes at zero energy. This is a direct 
consequence of the linear dispersion near K, in marked contrast to conventual two - 
dimensional electron gas, where the density of states is constant [78]. The unique band 
structure of graphene just described, brings about profound changes in the electronic 
properties of the system. In the continuum limit and in the effective - mass approximation, 
The Hamiltonian in the vicinity of the K point is [79] 


H.(K) a a oe ) a8 =-H5(K), (43) 


where @ = (ox, oy) is a vector of Pauli matrices. The Hamiltonian therefore describes 
two - dimensional massless neitrinos with a linear dispersion, ex(K) = tvrhk. The wave 


functions of these relativistic - like (Dirac) particles have a spinor structure, or the K and K’ 
points, they are without normalization 


aw, Ss i Pree er 
4 K = glk? — : _, & me Qik? 
Vi g*=e ( oe ) Wz =e ( : } (44) 


where s = +1 for upper band (electrons) and - 1 for the lowest band (holes0, tan6z = 
k,/k;. the upper and lower terms correspond to the quantum mechanical amplitudes or 
pseudospin) of finding the particle on one of the two sublattices A and B. In such 
approximation the Hamiltonian and wave functions for either valley can be written as the 2 x 
2 matrix 


— 0 k = ik =, és 
5 = : x y ’ wy 4 ‘3 = gik? : . 
H.(K) ome Kk, oo ft ( a ) (45) 


interestingly, in graphene, the pseudospln direction is associated with the momentum of 
the particles [80]. This means that the wave functions in the vicinity of the K and K’ points 
(Dirac points, Fig. 10) are chiral, or helical fermions. One consequence of this is that any 


backscattering, i.e scattering particles from the wave vector K to -K, is supressed (see, e.g. 
[86, 89]). Particles have opposite chirality in the K or K’ valley or in the electron or hole 
bands. The energy bands in the vicinity of these Dirac points are two cones meeting at k = 0 
(see also Fig. 6). The charge carrier in graphene are usually described as massless 
fermions [63, 79]. Experimental observation [77] of full Dirac cone dispersion have been 
measured in [81]. in such way the graphene provides unexpected connections between 
condensed matter physics and quantum field theory of the elementary particle physics. 
Phonons denote the quantized normal mode vibrations that strongly affect many 
processes in condensed matter systems, including thermal, transport and elastic properties 
(see, e.g. [91]). Phonons provide a sensitive probe of the electronic structure through 
Raman techniques since the absence of a band gap in graphene makes all incident 
wavelengths resonant, so that the Raman spectrum contains information about both atomic 
structure and electronic properties. An understanding of the phonon dispersion of graphene 
is essential to interpret the Raman spectra of graphene. Since the unit cell of monolayer 
graphene contains two carbon atoms, A and B, there six phonon dispersion bands (see Fig. 
8in chapter 5 in [15]), in which three are acoustic branches (A) and the other three are optic 
(O) branches. For one acoustic branch (A) and one optic (O) phonon branch, the atomic 
vibrations are perpendicular to the graphene plan, and they correspond to the out ~ of - 
plane(o) phonon modes. For two acoustic and two optic phonon branches, the vibration are 
in - plan (i). Traditionally, the directions of the vibrations are considered with respect to the 
direction of the nearest carbon - carbon atoms and, therefore the phonon modes are 
classified as longitudinal (L) or transverse (T) according to vibrations parallel with or 
perpendicular to, respectively, the A - B carbon - carbon directions. Therefore, along the 
high symmetry ['M and IK directions, the six phonon dispersion curves are assigned to LO, 
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iTO, oTO, LA, iTA, and oTA phonon modes. 
3. Particles phenomenology. 
3.1. Range of forces. 


The inability to deserve the stable of the atom by classical electrodynamics led to the 
creation of quantum electrodynamics (QED). 

A qualitative estimation of the radius of action of nuclear force can be made as follows. 
Since the nuclear forces do not depend on the electric charge and are of a purely quantum 
nature, the radius of action of the forces can only depend on mass of the carrier particles 
mM, and world constants # and c. Of these three quantities, only one constant of the 
dimension of length can be compiled - the Compton pion wave length 


R ~wZe, (46) 


Where m, is the pion mass If we assume that the pion travels at the speed of light ~ c, 
then in time z it travels a path 


R~cr~ bs. (47) 


The obtained distance is the radius of action of nuclear forces. 

Naturally, a complete quantum mechanical calculation of the radius of nuclear forces 
was made by Yukawa in 1935 [16]. 

The smallest entities of accessible strongly interacting matter in the world are hadrons, 
either baryons that are aggregates of three quarks or mesons that are made from quark - 
antiquark pairs. The most stable baryons, the protons and neutrons, are the major 
constituents of atomic nuclei, and the lightest meson is the pion. Understanding the 
structure of hadrons and how the properties of these particles arise from quantum 
chromodynamics(QCD) is a major interest of nuclear physics [83] (see, also [84]). This part 
of our review is devoted to pions and other strongly interacting bosons. For this reason we 
briefly describe the Yukawa approximation in task of strong interaction [85]. In 1935 Yukawa 
[18] tried to develop a theory of nuclear forces (see, e.g. [85]): The most important feature 
of these forces is that they have a small range. That is, nuclear forces decrease extremely 
rapidly when the interacting particles are more than about 10-5 m apart. The fundamental 
idea thus antedates the discovery of the pion by years. The role mesons in nuclear physics 
was not discovered experimentally; it was predicted through a brilliant theoretical 
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speculation. Below we shall introduce the Yukawa potential in its simplest form by analogy 
with the electromagnetic interaction. 

According quantum electrodynamics (QED) the scalar potential Ao produced by a 
charge distribution qp(x’) satisfies the wave equation [86]: 


ViAe- oe =-4nqp, (48) 


here q - electric charge and p is the nucleon probability density. If the charge distribution 
is time - independent, the wave equation reduces to the Poisson equation, 


V2Ao = ~ 4nqp. (49) 


in QED the simplest situation arises If the electromagnetic field is produced by a point 
charge q, at rest at x’. The potential is then given by 


Ao(X) = = (50) 


be 


and the interaction is ordinary Coulomb energy [87]. In this case, we have the solution 
of the Poisson equation in the next form: 


Ao(x) = [dx’ 22. (51) 


x-x'] ” 
For a point charge q located at the origin, Ao reduces to the Coulomb potential 
Ao(t) = +. (52) 


The important step at the solution of the Poisson equation can be summarized in the 
relation V2(+) = -426(x), where 6 is the Dirac delta function. 

Yukawa noticed that the electromagnetic interaction could provide a model but it did not 
fall off rapidly enough with distance. To force a more rapid decrease, he added a term k2 
to Eq. (49): 


(V? - K?)@(x) = 40 p(x). (53) 
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The last equation is the Klein-Fock-Gordon equation [86]. The electromagnetic potential 
Ao has been replaced by the field ®(x), and the strength of the field is determined by the 
strong source gp(x), where g determines the dimensionless strength, and p, as early, is 
probability density. The sign of the source term has been chosen opposite to the 
electromagnetic case. The solution of Eq. (53) that vanishes at infinity is 


is “k[x-X ’ ’ 
00) = arf BEE Drone. (64) 


For a strong point source, placed at position x’ = 0, this solution becomes the Yukawa 
potential [18], 


O01) = ar SS. (55) 


The constant k can be determined by considering Eq. (53) for the case (p(x) = 0) and 
comparing it to the corresponding quantized equation. For substitution 


E-ind, P=-inv, (56) 
changes the energy - momentum relation, 
E? = (pc)? + (mc?)’, (57) 


into Klein - Fock - Gordon equation 


Cay 
Sl 
Bho» 


2 Ve (42)? Joo =0. (58) 
For a time ~ independent field and p(x) = 0, comparison of Eqs. (57) and (58) yields 


k= 2. (59) 
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The constant k in the Yukawa potential is just the inverse of the Compton wavelength of 
the field quantum (see Table 2). The mass of the field quantum determines the range of 
the potential. The simple form of the Yukawa theory thus provides a description of the 
strong potential produced by a point nucleon in terms of the mass of the fleld quantum. It 
explains the short range of the strong forces. Below we briefly describe the most important 
properties of the nuclear force, that is, the strong force between a proton and an neutron. 


Table 4. Comparison between electromagnetic and Yukawa fields (after {88)). 


Electromagnetic field Nuclear field 


Static potential V=g/r V=g(e*/r) 
Eq. of propagation for potential v2V - ()(P Vat?) =0 VV-k2V- (4-\(e?Viat?) =0 
Relationbetween =f/panda ++ =0 24+ -k=0 
G CG 
. 2 2 f 
Relation between p=42 andE -+£>=0 t+ fy -k=0 
or E=cp (C?p? + A?k?¢2) 12 £ 


1. Attraction. The force is predominantly attractive, otherwise stable nuclei could not 
exist. 

2. Range and strength. The comparison of the binding energies of *H, 3H and “He 
(see, e.g. [88]) indicates that the range of the nuclear force is of the order 1 fm. If the force 
is represented by a potential with such width, a depth about 50 MeV is found [89]. 

3. Charge independence. The strong force is charge independent. After correction for 
the electromagnetic interaction (87, 90, 91] the pp, nn and np forces between nucleons in 
the same states are identical. 

4. Saturation. If every nucleon interact attractively with every other one, there would be 
A(A - 1)/2 distinct interacting pairs. The binding energy would be expected to be 
proportional to A(A - 1) = A?, and all nuclei would have a diameter equal to the range of the 
nuclear force. Both prediction, binding energy proportional to A? and constant nuclear 
volume, disagree violently with experiment for A ) 4. For most nuclei, the volume and the 
binding energy are proportional to the mass number [92]. Consequently, the nuclear force 
exhibits saturation: one particle attracts only a limited number others; additional nucleons 
are either not influenced or are repelled. A similar behavior occurs in chemical bonding and 
with van der Waals’ forces (see, e.g. [30]). Saturation can be explained through exchange 
forced or through strongly repulsive forces at short distances (hard core) [88]. 

5. Spin dependence. The force between two nucleons depends on the orientation of 
the nucleon spins [93 - 97]. 

6. Noncentral forces. Nuclear forces contain a noncentral component (see, also 
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below). 

To conclude this paragraph we should stress that the perturbation theory does not use 
for strong interaction. Furthermore, pions are not the only particles in the nucleon - nucleon 
force (Fig. 11). 
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Fig. 11. The diagrams corresponding to the fundamental pion - nucleon 
interactions. Charge independence requires that the coupling constants a, b, c, be 
connected a = (2b)!? and b = - c, b - The combining diagrams giving all cases of 
nucleon ~ nucleon scattering. The scattered amplitude for (p,n) is bc + a”, for (p,p) it 
is b?, for (n,n) it is c? (after [25]). 


As is well - known [88, 89] in 1947, the Yukawa particle, the pion was finally discovered 
in nuclear emissions. After 1947, more mesons kept turning up, and at present the list is 
long (see, e.g. [89] and references therein). 

As was noted above, once major problem in contemporary particle physics to explain 
why quarks and gluons are never seen as isolated particles. Usually, most particle 
physicists use "fields" and "particles" interchangeably, i.e. as denoting the same things. 
That is because the almost universal use of Feynman diagrams (see, e.g [98] and reference 
quoted therein) gives the false impression that particles (quanta) are always exchanged, 
even when they do not exist. The use of Feynman diagrams can be justified in weak 
nonlinear theories (weak coupling limitO but breaks down for strongly coupled nonabelian 
theories. This is correct also for strongly couple abelian theories with sources. This is an 
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example of mistaking the approximation (perturbation theory) for the exact theory. We try 
below to show that color confinement is a direct result of the nonabelian, i.e. nonlinear, 
nature of the color interaction QCD. This make it in general impossible to describe the color 
field as collection of elementary quants (e.g. gluons). In quantum field theory [99] an 
elementary particle, i.e. quantum = a harmonic excitation of a fundamental field (see, e.g. 
[86] , is defined through the creation and annihilation operators, a‘ and a, of the "second - 
quantized" theory. For instance, in QED, the entire electromagnetic field can be seen as a 
collection of superposed quanta, each with energy a:. the Hamiltonian of the 
electromagnetic field (omitting the zero - point energy) can be written 


H = ONkox (60) 
k 

where 

Nx = ala, (61) 


is the "number operator", i.e. giving the number of quanta with a specific four - 
momentum k when operating on a free state, 


Nic [eceefiees) = Mk [eee eMieee). (62) 


As all oscillators are independent, 


isoitikeee) 0 | Nx), (63) 


where nx is a positive integer, the number of quanta with that particular momentum. The 
energy in the electromagnetic field is thus the eigenvalue of the Hamiltonian (64). The 
reasoning for fermion fields is the same, but then the number of quanta in any given state 
can be only 0 or 1 (Fermi statistics). 

Now, assuming that QCD is the correct theory of quark interactiens, a problem arises, 
as it is general impossible to writeof the cclor fields in termsf superimposed harmonic 
oscillators. It is not possible to represent the solution as a Fourier expansion and then 
interpret the Fourier coefficient [100] as creation/annihilation operators through "second 
quantization", as the color vector potentials A® (b < 1, ...,8) are governed by nonlinear 
evolution by nonlinear evolution equations, 
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DY Fw = jy, (64) 


and Fourier methods are inapplicable to nonlinear equations [100]. 
Without a quark current, j, = Gs¥y,yw = 0, we get, in component for 


(Sab0" + QsfarcAt)(,A? - OVA + OsfodeAZAy) = 0, (65) 


where gs is the color coupling constant (Summation over repeated indices implied). 
When we have an abelian dynamical group, as in QED, all the structure constants fap, are 
zero. Eq. (65) then reduces to a linear differential equation, and a general solution can be 
obtained the Fourier expansion 


AGE? = fae 5 [ax(A)en(k, 2)e™ + al (Aes (k, 2) eq, (66) 
A=0 


where g, is the polarization factor. 
However, for a theory based on a nonabelian group, like QCD, this can no longer be 
done, due to the nonlinear nature of Eq. (65) when fabcfrde # 0, 


3 
ARO + Fok D7 [aP(A)en(k, Ae + abi(ayen (k, Ae", (67) 
A =O 


Thus, the color fields can be represented by harmonic oscillators (gluons) only in the 
trivial, and physically empty, limit when the strong interaction coupling constants tends to 
zero, gs — O (or, within perturbation theory, equivalently when Q? — oo because of 
asymptotic freedom. Hence, no elementary quanta of the color interaction, in the usual 
sense, can exist. This means that no gluon particles are possible, and that Eq. (60) does 
not hold for color fields. The fields are there, but their quanta, gluons (and quarks), are 
relevant only when probed at sufficiently (infinitely) short distances. we note that a quark 
field is the source of a color field, but this color is itself a source of a color fie;d. Hence, a 
quark field is never removed from other sources, is always interacting, and can never be 
considered _to be freely propagating. This results in that the quark fields can never be 
represented by harmonic oscillator modes, unless Q? -+ oo, whereas physical particles, 
observable in nature, should exist as Q? — o. this means that no quark field quanta 
(quarks) can ever exist. In QCD, the particles "giuons" and "quarks" are merely artifacts of 
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the approximation method used, i.e. the perturbative expansion in the interaction on a 
"background" of assumed free gluons and quarks. They are simply absent in the exact 
theory (in details see [101]). 

Contemporary theory predicts that both quarks and gluon acquire running mass 
distribution in QCD, which are large at infrared momenta [15]. The current - quark of 
perturbative QCD evolves into constituent - quark as its momentum becomes smaller. The 
constituent - quark mass arisés from a cloud of low - momentum gluons attaching 
themselves to the current quark. This is dynamically chiral symmetry breaking (DCSB) : an 
essentially nonperturbative effect that generates a quark mass from nothing, namely it 
occurs even in the chiral limit. DCSB namely the generation of mass from nothing, is a 
theoretically - established nonperturbative feature of QCD [14]. 


3.2. The quark ~ gluon interaction. 


lt is well - known that including the color property, a kind of quark charge, the Pauli 
principle may be asserting. As was shown above, the quantum number color can assume 
three values, which called red, blue and green. Accordingly, antiquark carry the anticolors 
anti - red, anti - blue, and anti - green. Now the three u - quarks (for example in jA**) = 
lu f u fu f)) may be distinguished [102]. Thus, a color wave function antisymmetric under 
particle interchange can be constructed, and we so have antisymmetry for the total wave 
function, We all know too that the interaction binding quarks into hadrons is called the 
strong interaction. Such a fundamental interactions always connected with a particle 
exchange. For the strong interaction, gluons are the exchange particles that couple to the 
color charge of quarks. This is analogous to the electromagnetic interaction in which 
photons are exchanged between electrically charged particles. Tne modern experimental 
findings (see, e.g. [103]) led to the development of a field theory called quantum 
chromodynamics - short QCD. As its name implies, QCD is modelled upon quantum 
electrodynamics - QED [104]. In both, the interaction is mediated by exchange of a 
massless field particle with JP = 1° (a vector boson). Each gluon carries one unit of color 
and one of anticolor. It should appear, then, that there should be nine species of gluons - rr, 
rb, rg, br, bb, bg, gr, gb, gg. Such nine - gluon theory is perfectly possible in principle, but it 
would describe a world very different from our own [105]. In terms of group theory of color 
SU(3) symmetry (on which QCD is based [14, 106]), the 3 x 3 color combinations form two 
multiplets of states: a ’coloroctet’: 
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[ 11) = (rb + br)//2 15) = -i(rg - gf)//2 

I2) = ~i(rB - bt) 2 16) = (b9 + gb) /2 
13) =(rF- bb)//2 I7) =-i(b9 - gb)/ 2 
4) = (79+ gf\//2 18) =: (rf + BD - 298)/,/6 


(68) 


and a ‘color singlet’ 


19) = (rf + bb + 98)1//3. (69) 


If the singlet gluon existed, it would be a common and conspicuous as photon. May be 
the ninth gluon is the photon. Confinement requires that all naturally occurring particles be 
color singlets, and this explain why the octet gluons never appear as free particles [107, 
108]. But 19) is a color singlet, and if it exists as a mediator it should also occur as a free 
particle. Moreover, it could be exchanged between two color singlets (a proton and a 
neutron, say) giving rise to a long - range force strong coupling, whereas in fact up to 
present time we know that the strong force is of very short range. Therefore, there are 
evidently only eight kinds of gluons. We should add because gluons are massless, they 
mediate a force of infinite range ( the same as electrodynamics). In this since the force 
between two quarks is actually long range. However, confinement, and the absence of a 
singlet gluon, conceals this from us. A singlet state (such as the neutron) can only emit and 
absorb a singlet (such as the pion), so individual gluons cannot be exchanged between a 
proton and neutron. That’s why the force we observe is of short range. If the singlet gluon 
existed, it could be exchanged between singlets, and the strong force would have a 
component of infinite range. 

By their exchange the eight gluons mediate the interaction between particles carrying 
color charge (quarks), i.e. not only quarks but also the gluons themselves. This is an 
important difference to the electromagnetic interaction, where the photon field quanta have 
no charge, and therefore cannot couple with each other. In analogy processes of QED 
(emission and absorption, pair production and annihilation); emission and absorption of 
gluons (Fig. 122 take place in QCD, as do production and annihilation of quark - antlquark 
pairs (Fig. 12°). In addition, however three or four gluons can couple to each other in QCD 
(Fig. 12°), 
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a) b) 3 ) qd) 


Fig. 12. The fundamental interaction diagrams of the strong interaction: emission 
of a gluon by a quark (a), splitting of a gluon into quark - antiquark pair (b), and ’self - 
coupling’ of gluons (c, d). 


Further we briefly consider some main properties of group SU(3). The set of unitary 3x3 
matrices with det U = 1 form the group SU(3). The generators may be taken to be any (3? - 
1 = 8) linearly independent traceless Hermitian 3x3 matrices (see e.g. [14]). Since it is 
possible to have only two of these traceless matrices diagonal, this is the maximum 
number of mutually commuting generators. This number is called the rank of the group 
[106], so that SU(3) has rank 2, and SU(2) has rank 1. According to the theory of groups, 
the fundamental representation of SU(3) is a triple. The three color charges of a quark R, G, 
and B form the fundamental representation. of an SU(3) symmetry group. in this 
representation, the generators are 3x3 matrices, which is called by A;, where | = 1, 2 ....8. 
For Ai, the convention is to adopt Gell-Mann matrices [109], which are defined as 


010 0 -i 0 100 001 

At = 100 A2 = ioo Ag={ 0-10 Ag = 000 
000 000 000 100 
00 -i 000 000 100 

As=] 000 4e=| 001 ar=] 00-4 dg = 010 
i00 010 010 00 -2 


(70) 


with simultaneous eigenvectors 
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The eight 4;s satisfy the following commutations relations: 


A 4 
[ 4, a = ify (72°) 
[2, Z] = Foy + dye, (73°) 


where fj is called the structure constant of SU(3) with {a, b} = ab + ba and fy, dix are 
totally antisymmetric or symmetric with respect to i i k permutations and have the values 
listed in Table 5. 


Table 5. Structure constant of SU(3). 


frog = 1 


1 
f147 = fags =fos7 = fags = fs16 = fea7 = > 


ED 
fase =fe7e = -5- 


= = a = LL 
diis = d22g = dag = - dees 5 


dias = d157 = dase = daa4 = dass = + 


dea7 = ses = 377 = - + 


d = d = d = d = aa ee 
448 558 668 779 aie 


Below we briefly consider the interaction of a quark and anti - quark, in QCD. We shall 
assume that they have different flavors, so the only diagram (in lowest order is one in Fig. 
13)., representing for instance , us d -- u+d. The amplitude of this process is given (see 
Fig. 13) by (see e.g. [110] 


M = “3 S[u(8)y"u(t)JIv(2)y ww(4)](0§4%01 )(C2A%C4) (74) 


Here summation over @ implied. This is exactly what we have for electron - positron 
scattering in QED, except that we use g; (of course), and we have in addition the ‘color’ 
factor 
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= L(ChA%c1 )CsA%C4) . (75) 


The potentlal describing the qq interaction is, therefore, the same as that acting in 
electrodynamics between two opposite charges, only with « replaced by fas: 


Vag(r) = - F922. (76) 


Now, the color factor itself depends on the color state of the interacting quarks. From a 
quark and an antiquark we can make a color octet (68) and a color singlet (69) (all members 
of which yield the same f). A typical octet state (68) is rb (any of the others would do just as 
well) .Here the incoming quark is red, and outgoing quark must also be red and the 
antlquark antiblue. Thus 


1 0 
C1 =C3= QO |, Co=C4= 1 (77) 
4) 0 
1 0 
and hence f= [(100)a"} 0 [H(010)A*] 1) |] = 4A4,A5p. (78) 
0 0 
3,23 ; Pq, &q 
> 
Py,0, 


Px, 0g 


Fig. 13. The quark - antiquark interaction. 


A glance at the 2 martices (70) reveals that the only ones with entries in the 11 and 22 
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positions are A3 and Ag. So we can write 
f= L(Ag Pap?) + al ae) = 211-1) + (1S) (1/73 ) 1 =- $. (79) 


The color singlet state is (69) [(e# + bB + g)1//3] If the incoming quarks are in the 
singlet state the color factor is a sum of three terms: 


1 0 0 

ses | 447] 0 | [(100)Asca+| cae] 1 (010)A%c4+| e442] 0 | |(001)A%¢ 
0 0 1 

(80) 


The outgoing quarks are necessary also in the singlet state, and we get nine terms in all, 
which can be written compactly as follows: 


f= $5 AR (ayy) = THA) (81) 


in last formula summation over i and j, from 1 to 3, implied in the second expression. 
Now 


Tr(AZA") = 26%, (82) 
so, with the summation over «, 
THAA") = 16 (83) 
Evidently, then, for the color singlet 
f=4 (84) 


Putting equations (80) and (85) into equation (77), we conclude that quark - antiquark 
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potentials are 
Van(t) =- 4 {42 (85) (color singlet) and 
Vaa(t) = $2. (86) (color octet). 


from the signs we see that the force is attractive in the color singlet but repulsive for the 
octet. This helps to explain why quark - antiquark binding (to form mesons) occurs in the 
singlet configuration but not in the octet (which would have produced colored mesons) (see, 
also [11]) . 


3.3, Modern view of the origin of the nuciear force. 


As is well - known the nuclear force is that force which binds together the nucleons in 
the nucleus. Historically, the nuclear force was called the Strong Nuclear Force and was 
considered to be one of the four forces, along with the gravitational force, the 
electromagnetic force and the weak nuclear force [111, 112}. Soon after the discovery of 
quarks, the force which holds the quarks together in a nucleon was called the Cromo 
(hadron) Dynamic Force. This chromo dynamic force, which is a sub ~ nucleon force, is 
considered to be much larger than the nuclear force. Later, the chromo dynamic force was 
redefined as the strong nuclear force, which is now considered to be a sub - nucleon force 
responsible for the behavior and interactions of sub - nucleon particles. The force which 
hold the nucleons together in a nucleus was renamed as the Nuclear Force. to add to this 
confusion even more, there is a model of the nuclear force, which is called the residual 
chromo dynamic force. Because of this model, it is presumed that the nuclear force is 
simply a subset of the strong nuclear force. Because of this assumption, it is still claimed 
that there are only four forces in nature: thestrong nuclear force, the electromagnetic force, 
the gravitational force, and the weak nuclear force. So The strong nuclear force has two 
parts, the chromo dynamic force which is sub - nucleon force, and the nuclear force, which 
is that force holding the nucleons together in a nucleus. tn this part of our chapter we use 
the term nuclear force, which will be used to describe that force which holds together the 
nucleons in a nucleus. Our present view is that nuclear force may be divided into three 
parts, as illustrated schematically in Fig. 14. The long - range part (r » 2 fm) is dominated by 
one - pion exchange. !f exchanges of a single pion are important, there is no reason to 
exclude similar processes involving two or more pions and mesons heavier than pions. 
The range of interaction associated with these more massive bosons is shorter, and for this 
reason, the intermediate - range part of the nuclear force (1fm < r <« 2 fm) comes mainly 
from exchanges of single heavier mesons and two pions. The hard core in the interaction (r 
<1 fm) is made of heavy meson exchanges, multipion exchanges, as well as QCD effects 
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(the more detail see [3]) 


Yo 


‘hard scalar meson Bien 
core exchange exchange 


Fig. 14, Schematic diagram showing the different parts of a nucleon - nucleon 
potential as a function of distance r between two nucleon. The hard core radius is 
around 0.4 fm and it takes energy > 1 GeV to bring two nucleons closer than (twice) 
this distance. The main part of the attractions lies at intermediate ranges, at radius ~ 
1 fm, and is believed to be dominated by the exchange of scalar mesons. The long - 
range part, starting at around 2 fm, is due to the single - pion exchange. 


Now we turn to the task of understanding the strength and the form of the nuclear force 
from the structure of the nucleons and the strong interaction of the quarks inside the 
nucleons. Naturally in the following discussion we will employ qualitative arguments. The 
structure of the nucleon will be approached via the nonrelativistic quark model where the 
nucleons are built out of three constituent quarks. The nuclear force is primarily transmitted 
by quark - antiquark pairs, which we can only introduce ad hoc through plausibility 
arguments. We should repeat that a consistent theory of the nuclear force, based upon the 
interaction of quarks and gluons, does not yet exist [10, 88, 92, 112]. 

We begin with the short distance repulsive part of the nuclear force (see Fig. 14) and try 
to construct analogies to better understood phenomena. We remember that atoms repel 
each other at short distances is a consequence of the Pauli principle. The electron clouds 
of both atoms occupy the lowest possible energy levels and if the clouds overlap then some 
electrons must be elevated into excited states using the kinetic energy of the colliding 
atoms. Hence we observe a repulsive force at short distances (see, also [78]). The quarks 
in a system of two nucleons also obey the Pauli principle, i.e. the 6 quark wave function 
must be totally antisymmetric. It is, however, possible to put as many as 12 quarks into the 
lowest / = O state without violating the Pauli principle, since the quarks come in three colors 
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and have two possible spin ({, t) and isospin (u - quark, d - quark) directions. The spin - 
isospin part of the complete wave functions must be symmetric since the color part is 
antisymmetric and, for / = 0, the spatial part is symmetric. We thus see that the Pauli 
principle does not limit the occupation of the lowest quark energy levels in the spatial wave 
function, and so the fundamental reason for the repulsive core must be sought elsewhere. 
The real reason is the spin - spin interaction between the quarks [94, 95 113, 114]. The 
potential energy then increases if two nucleons overlap (see below Fig. 10) and all 6 quarks 
remain in the /= 0 state since the number of quark pairs with parallel spins is greater for 
separated nucleons. Of course the nucleon - nucleon system tries to minimize its 
*chromomagnetic’ energy of maximizing the number of antiparallel quark spin pairs. But this 
is incompatible with remaining in an /= 0 state since the spin - flavor part of the wave 
function must be completely symmetric. The colormagnetic energy can be reduced if at 
least two quarks are put into the /= 7 state. The necessary excitation energy is comparable 
to the decrease in the chromomagnetic energy, so the total energy will in any case increase 
if the nucleons strongly overlap. Hence the effective repulslon at short distances is in equal 
parts a consequence of an increase in the chromomagnetic and the excitation energies. If 
the nucleons approach each other very close (r = 0) one finds in a non - adiabatic 
approximation that there is an 8/9 probability of two of the quarks being in a p state (see, 
e.g. [115, 116]. This configuration express itself in the relative wave function of the nucleons 
through a node at 0.4 fm. This together with the chromomagnetic energy causes a strong, 
short range repulsion. The nuclear force may be described by a nucleon - nucleon potential 
which rises sharply at separation less than 0.8 nm [89, 92]. 

Further we turn to the attractive part of the nuclear force. Again we will pursue analogies 
from atomic physics. As we know the bonds between atoms are connected to a change in 
their internal structure and we expect something similar from the nucleons bound in the 
nucleus. Indeed a change in the quark structure of bound nucleons compared to that of 
their free brethren has been observed in deep inelastic scattering off nuclei (EMC effect 
(117, 118]). It is clear upon a moments reflection that the nuclear force is not going to be 
well described by an ionic bond [78]: the confining forces are so strong that it is not possible 
to lend a quark from one nucleon to another. A van der Waals force, where the atoms 
polarize each other and then stick to each other via the resulting dipole - dipole interaction 
can also not serve us as a paradigm (see, also [119, 120]). A van der Waals force 
transmitted by the exchange of two gluons (in analogy to two photon exchange in the atomic 
case) would be too weak to explain the nuclear force at distances where the nucleons 
overlap and confinement does not forbid qluon exchange. At greater separations gluons 
cannot exchanged because confinement. Although color neutral gluonic states (glueballs) 
could still be exchanged, none which are light enough have ever been experimentally 
observed (see however below experimental manifestation of the strong interaction in solids). 

The only analogy left to us to explain the nuclear force is a covalent bond [78], such as 
that which is, e.g., responsible for holding the Hz molecule together [121]. Here the 
electrons of the two H atoms are continually swapped around and can be ascribe to both 
atoms. The attractive part of the nuclear force is strongest at distances around 1 fm and 
indeed reminds us of the atomic covalent bond. To simplify what follows, let us assume that 


~ 406 - 


the nucleon is made up of a two quark system (diquark”) and a quark (Fig. 15). Such 
description has proven to be very successful in describing many phenomena [16]. 


Nucleon Nucleon 


Fig. 15. Quark configurations in a covalent bond picture. At large separations, 
when the nucleons just overlap, we may understand them as each being diquark - 
quark systems. 


The most energetically favorable configuration is that where a u - and d - quark combine 
to form a diquark with spin 0 and isospin 0. The alternative spin 1 and isospin 1 diquark is 
not favored. The covalent bond is then expressed by the exchange of the single quarks. 
Since the nuclear attraction is strongest at distances of the order of 1 fm [9, 85], we do not 
to worry about confinement effects. The covalent bond contribution to this force can be 
worked out analogously to the molecular case. However, the depth of the potential that is 
found in this way is only about one third of the experimental value [122, 9, 33]. In fact quark 
exchange is less effective to be exchanged the quarks must have the same color, and there 
is only a 1/3 probability of this. The contribution of direct quark exchange sinks still further if 
one takes the part of the nucleon wave function into account where the diquarks have spin 1 
and isospin 1. Thus the covalent bond concept, if it is directly transferred from molecules to 
nuclei, does not give us a good quantitative description of what is going on in nuclei. It 
should be noted that this is not a consequence of confinement, but rather of direct quark 
exchange being suppressed as a result of the quarks having three different color charges. 
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*) See, also: S.-X. Qin, C.D. Roberts, S.M. Schmidt, Spectrum of light - and heavy - 
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baryons, Few - Body Systems 60, 26 - 44 (2019); M.Y. Baranov, M.A. Bedola, W.K. Brooks 
et al., Diquark correlation in hadron physics: Origin, Impact and Evidence, ArXiv/ hep - ph/ 
2008.07630 (2020). 


One more the degree of freedom is an effective quark - quark exchange which may be 
produced by color neutral quark - antiquark pairs (see quarks). These pairs are continually 
being created from gluons and annihilated back into them again. This quark - antiquark 
exchange actually plays a larger role in the nucleon ~ nucleon interaction than does the 
simple swapping of two quarks. It must be stressed that this exchange of color neutral quark 
- antiquark pairs does not only dominate at great separations where confinement only 
allows the exchange of color neutral objects but also at relatively short distances. One may 
thus understand the nuclear force as a relativistic generalization of the covalent strong force 
via which the nucleons finally exchange quarks (the more details see [9, 10, 16}). 


3.4. The Concepts of Quantum Chromodynamics. 


Quantum chromodynamics, familiarly called QCD, is the modern theory of the strong 
interaction (see, e.g. [11, 16]). Historically its roots are in nuclear physics and the 
description of ordinary matter - understanding what protons and neutrons are‘and how they 
interact. QCD is renormalizeable quantum field theory of the strong interaction. Its 
fundamental constituents are spin 1/2 Fermi fields called quarks carry fractional electric 
charge and non - abelian spin gauge fields called gluons which interact with the quarks as 
well as among themselves. QCD is the realization that simple principle of a local gauge 
group SU (8) and the formal existence of fermions called quarks transforming as the 
fundamental representation of this group essentially specifies the complete theory. The 
main idea of QCD is to make the SU (3) color symmetry a local, rather than just a global 
symmetry [123]. As is well known, to implement this local SU (3) symmetry, one must 
introduce a compensating vector gauge field [124] A#(x), a= 1, 2,... 8 which transforms 
as an adjoint representation of SU (3). The octet of gauge fields are called gluons (they do 
not carry flavor). 

Thus, the theory that describes strong interaction in the Standard Model [125] is called 
QCD. Although this is not tested to the same extent or precision as QED, it is nevertheless 
in impressive agreement with large body of experimental data and is not contradicted by 
any known experiment. QCD is similar QED in that both describe interactions that are 
mediated by massless spin 1 bosons coupling to conserved charges. Theories of this type 
are called gauge theories because they have a characteristic symmetry called gauge 
invariance” 


Ane Sea Ae mae Sm nd Melt eA er UY es SR? crm init Ni 


*) In the gauge principle we use [126] 
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E=v¢-2, Be¥xK 1) 

defined by ¢(7, t) and A(?, t) are the scalar and vector potentials. In particular, one 
easily shows that substititing 

dog =64+4, RiW-=R-F EQ) 

in (1’) leaves the electric and magnetic fields unchanged, when f(7, t) is an arbitrary 
scalar function. The replacement of (¢, A) by (0. x) in accordance with (2’) is called a 


gauge transformation, and a theory whose physical predictions are unaltered by such a 
transformation is said to be gauge - invariant. Because observable quantities depend only 
on the fields E and B, and not on the particular choice of potentials used to describe them, it 
is fundamental requirement that any theory formulated in terms of these potentials should 
be gauge - invariant. 


As is well known the gauge invariance plays a fundamental role in theoretical treatments 
of QED and QCD, where it can be used to infer the detailed forms of the interactions. Here 
we shall adopt a more phenomenological approach. The spin 1 bosons are called gauge 
bosons. In QED they are photons, and in QCD they are called gluons. Gluons have zero 
electric charge, like photons, but couple to the color charges that rather than to the electric 
charge. This leads immediately to the so - called flavor independence of strong interactions: 
i.e. the different quark flavors a = u, d, s, c, b and t must have identical strong interactions, 
because they exist in the same three color states r, g, b with the same possible values of 
the color charges. This has its most striking consequences for u and d quarks, which have 
almost equal masses, where it leads to the phenomenon of isospin symmetry [14]. 

A second property of strong interactions that follows from the above picture without 
detailed argument is that the forces between the quarks must be long range, because the 
gluons have zero mass. This does not imply that the forces between hadrons are also long, 
because hadrons have zero color charges overall. The forces between the colorless 
hadrons are the residues of the forces between their quark constituents, and cancel when 
the hadrons are far apart”. 
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*) The main question are: How long - range electromagnetic quarks interaction inside 
hadrons reduces on the borders of hadrons to short - range nucleons strong interaction. 


We have noted that QED and QCD both describe interactions, albeit of very different 
strengths, which are mediated by massless spin 1 bosons that couple to conserved 
charges. However, there is a crucial difference between them which profoundly affects the 
character of the resulting forces. It is that while the photons that couple to the electric 
charge are themselves neutral, gluons have nonzero values of the color charges to which 
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they couple. This is illustrated in Fig. 16,whilch shows a particular example of a quark - 


quark interaction by gluon exchange, where the gluon is represented, as early, by a 
‘corkscrew’ line to distinguish it from a photon. 


u,b 
ur 


sb oF 


Fig. 16 Example of quark - quark scattering by gluon exchange. In this diagram 


the quark flavor u or s unchanged on gluon emission, but the color state can change, 
as shown. 


In this picture, the color states of the two quarks are interchanged, and the gluon has the 
color quantum numbers 


IF = 18(r) - 1$(b) = F (87) 
and 
Yo = ¥C(r) - Yb) = 1, (88) 


which follow from applying color conservation at, for example, the lower vertex and using 
the values for the quark color states r, b listed in Table 4. 


Table 4. Values of color isospin charge If and color hypercharge Y© for the color 
states of quarks and antiquarks. 
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Quarks Antlquarks 
ig ye aig Ye 
1/2 18 if -1/2 -1/3 
g -1/2 1/3 yg 1/2 -1/3 
b 0 -2/8 b- 0 28 


Just as quarks exist in three different color states, gluons can exist in eight different 
color states, which is determined by symmetry [11, 124], and the If, Y° values correspond 
to just one of them [15]. Further, we note that is gluons couple to particles with nonzero 
color charges, and if gluons themselves also have nonzero charges, then by implication 
gluons couple to other gluons. The two types of self - coupling that occur in QCD are 
illustrated in Fig. 17. The first is a gluon exchange process analogous to Fig. 17 for quark - 
quark scattering, while the second involves a zero - range ‘contact’ interaction (Fig. 17°). 
The gluon - gluon interaction of Fig. 17 have no analogue in QED, and it can be shown that 
they lead to properties of the strong interaction that differ markedly 


é 
@) (b) 


Fig. 17. The two lowest - order contribution to gluon - gluon scattering in QCD. 


from those of the electromagnetic interaction. These properties are color confinement 
and asymptotic freedom. Color confinement is the requirement that observed states have 
zero color charges. This was discussed early, where, for example, it was shown that all 
quark bound states (hadrons) must have integral electric charges. It also implies that 
gluons, like quarks, cannot be observed as isolated free particles, because they have 
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nonzero values of the color charges. Bound states of two or more gluons with Zero color 
charges overall can be formed in principle, due to the strong interaction between gluons 
themselves. such states are called glueballs and will be briefly discussed below. 

Asymptotic freedom means that the interaction get weaker at short distances, and at 
distances less than about 0.1 fm the lowest - order diagrams dominate. At these distances 
quark - quark scattering, for example, is given approximately by one - gluon exchange 
diagrams like Fig. 17. However, as the distance between the quarks increases, the 
interaction gets stronger and many high - order diagrams become important. In this strong 
interaction regime the situation is very complicated, and it has not yet been possible to 
evaluate the theory precisely. We therefore have to rely an approximate results obtained by 
numerical simulations of the theory on very large computers, and the demonstration of 
confinement in QCD rests largely on such simulations. As will be shown below, another 
important consequence from asymptotic freedom is the fact that the strong coupling as is 
small enough, at sufficiently large momentum transfers, to allow application of perturbation 
theory In order to provide quantitative predictions of physical processes (see, also [102, 
126)). 

The above features are conveniently illustrated by considering the static potential 
between a heavy quark and its antiquark in a color singlet state. This is appropriate potential 
for a nonrelativistic discussion of charmonium and bottomium, and was determined 
empirically for the limited range 0.2 < r < 0.8 fm [89]. Here we are concerned with the 
behavior predicted on the basis of QCD. At short distances r < 0.1 fm, the interaction is 
dominated by one - gluon exchange and we might expect a Coulomb - like potential 
analogous to that arising from one - photon exchange in QED. In fact it can be shown that 
the potential is given by 


Vir)=- 4% = (r < 0.1fm), (89) 


where the strong coupling constant as is a measure of the strength of the interaction 
analogous to the fine structure constant a in QED. Because of asymptotic freedom, the 
strength of the interaction, and hence a;, decreases with decreasing r, but for r < 0.1 fm this 
variation is slight and can in many application be neglected. At distances beyond 0.1 fm, 
however, the strength of the interaction increases more rapidly and one - gluon exchange 
no longer dominates. In this region, we are forced to rely on lattice gauge theory 
calculations of limited precision ( see, also [4]). These are consistent with the empirical 
potential in Fig. 18 in the range where this is determined, and show that at large distances 
the potential rises approximately linearly 


Virye ar (r > 1fm), (90) 


where the constant 4 cannot be calculated precisely, but is of order 1 GeVfm~!. This is 
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an example of confining potential in that it does not die away with increasing separation and 
the force between quark and antiquark cannot be neglected, even when they are very far 
apart. 
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Fig. 18. Heavy quark - antiquark potentials obtained from fitting the energy levels 
of charmonium and bottomium. The solid line show the results of calculation used 
next relation V(r) = -2 + br ( with a = 0.30 and b = 0.23 GeV?) and dashed line 
obtained used the relation V(r) = aln(br) (with a = 0.75 GeV and b = 0.80 GeV). 


To conclude this part we should highlight that QCD, the gauge field theory of the strong 
interaction, has specific features, asymptotic freedom and confinement, which determine 
the behavior of quark and gluons in particle reaction at high and at low energy scales. 


3.4.1. The strong coupling constant. 


In this part of our review, we give more formal discussion of the strong coupling ag in the 
perturbation domain. As is well - known the strong interaction, binding quarks and gluons 
inside hadrons, is the strongest of the four fundamental forces in Nature which we know 
today. In opposition to gravitational and electromagnetic forces, the strong force, as well as 
weak interaction, only acts at subatomic distances (see, however below). The strong force 
not only determine the binding of quarks and gluons inside hadrons, it also determines the 
cohesion of protons and neutrons inside atomic nuclei [14, 89]. Hadrons like protons and 
neutrons are responsible for more than 99 % of the mass of all visible matter in our 
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Universe, and those masses are mainly generated by the strong binding of quarks inside 
hadrons, rather than by the (generally small ) masses of the quarks themselves. It is a 
common place in contemporary physics that the restriction of the strong force to subatomic 
distances is a consequences of two characteristic features: they are called "Confinement" 
and "Asymptotic Freedom". Confinement is a necessary requirement to explain the fact 
that no isolated quarks have ever been observed in any experiment, although symmetry 
arguments and scattering experiments in the 1960’s established quarks [9, 84] with -1/3 or 
+2/3 of electrical charge units and a newly introduced quantum property called "colored 
charge", as the basic constituents of hadrons. Confinement determines that at large 
distances, or equivalent - at low momentum or energy transfer in elementary particle 
reactions. The strong force prevents the existence of free quarks: trying to separate two 
quarks from each other, for instance in high energy scattering reactions, apparently results 
in an increase of the force field’s energy at large distances, such the new quarks are 
created out of the vacuum - the initial quarks "dress" up with other quarks to build hadrons. 
These hadrons exhibit no net color charge to the outside, such that they appear as 
elementary entities rather than the quarks themselveseww3. 

The term "Asymptoti Freedom" is used to describe the behavior to describe of quarks at 
high energy or momentum transfer, or - equivalently - at small distances. Also this feature is 
based on experimental observations: in high energy scattering processes between leptons 
(e.g. electrons or neltrinos) with protons or neutrons, the dynamics reveal that scattering 
occurs at pointlike and massless constituents, the quarks, rather than at a homogeneous 
object with size of proton (neutron). Apparently, at sufficiently high momentum transfers, 
quarks behave like free of weakly bound particles. Also, quarks knocked out of a hadron, in 
a high energy scattering process, were never observed as free particles. this picture gives 
physical insight into the phenomenon of asymptotic freedom ("ultraviolet" UV regime) at 
high momentum transfer (at small distances) and to s strong coupling regime in the low 
momentum transfer regime ("infrared" IR regime ) (large distanes). Sometimes the 
confinement, also called "infrared slavery" [127 - 130]. 

The coupling strength as is the basic free parameter of QCD, the theory of strong 
interaction [16]. QCD describes the interaction of quarks through the exchange of an octet 
of massles vector gauge bosons, the gluons, using similar concepts as known from QED. 
While QCD does not predict the actual size of as at a particular energy scale, its energy 
dependence is precisely determined. If the renormalized coupling as(u?) can be fixed (i.e. 
measured) at a given scale »?, QCD definitely predicts the size of a, at any other energy 
scale Q? through the renormalization group equation [63] 


Q2 2. = plas(Q?)) (91) 


The perturbative expansion of the # function, including higher order loop correction to 
the bare vertices of the theory, is calculated to complete 4 - loop approximation [126]: 
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where N, is the number of active quark flavors at the energy scale Q. The numerical 
constants in equation (93) are functions of the group constants C,= N and Cr = 
(N? - 1)/2N for theories exhibiting SU(N) symmetry [55, 61]; for QCD and SU(3), Ca = 3 
and Cr = 4/3 (for details see [4, 55]. 

A solution of equation (92) in 1 - loop approximation, i.e. neglecting 81 and higher order 
terms, is 


re cs (94) 
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Apart from giving a relation between the values of as at two different energy scales Q? 
and 47, equation (94) also demonstrates the property of asymptotic freedom: if Q2 becomes 
large and Bo is positive , i.e, if Ny « 17, as(Q?) will asymptotically decrease to zero. 
Likewise, equation (94 ) indicates that as(Q?) grows to large values and, in this 
perturbative form actually diverges to infinity at small Q?: for instance, with at (2? = Mz, ) 


= 0.12 and for typical values of Ny = 2 ... 5, as(Q?) exceeds unity for Q? « 
6(100 MeV .... 1 GeV) Clearly, this is the region where perturbative expansions in as are 
not meaningful anymore, and we may regardenergy scales below the order of 1 GeV as the 
nonperturbative region where confinement sets in, and where equations (93) and (94) 
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cannot be applied. 
Including £; and higher order terms, similar but more complicated relations for as(Q?), 
as a function of as(2) and of Ing as in equation (94) emerge. They can solved 


numerically, such that for a given value of as(u?), choosing a suitable reference scale like 


the mass of the Za bozon, yw = Mz,, as(Q?) can be accurately determined at any energy 
scale Q? > 1 GeV?. 
If we set 


Qo oa 
e(Pors(x?)) ? 


a dimensional parameter A is introduced such that equation (94) transforms into 


as(Q?) = Fan CO (95) 


Hence, the A parameter is technically identical to the energy scale Q where a;(Q?) 
diverges to infinity as(Q?) — oo for Q? —+ A?. Further we give a numerical example [126] A 
~ 0.1 GeV for as(Mz, = 91.2 GeV) = 0.12 and Ny =5. 


In complete 4 - loop approximation and using the A - parametrization, the running 
coupling is thus given by (126, 128] by 


as(Q*) = aor ~ ar Asink + 

ate (2p intL - int - 1 ]+#) + 

BBLS \ AB Be 

—t-{ 2. in8b 4 Stn? ~ 1). 3 ibe, Bs 
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where L = Q?/Agg. The first line of equation (118) includes 1 - and the 2 - loop 
coefficients, the second line of the 3 - loop and the third line is the 4 - loop correction 
respectively. 

Measured values of as(s) [127 - 130], obtained from a variety of different processes are 


shown in Fig. 19, where the curves show the predicted behavior corresponding to the ’best - 
fit’ value [129] 
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as (Mz) = 0.118 + 0.002 (97) 
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Fig. 19. Values of the running coupling constant a; obtained in different 
measurements. 


at the reference value so = M, (see, also [128]). The decrease in as() as » increases, 
corresponding to shorter distances, is clearly seen, and because of this variation, as(s) is 
often referred to as the running coupling constant” (see, also [127 - 130]). 

With QCD and the electroweak theory in hand, what remains to be understood? If both 
theories are correct, can they also be complete? Many observation are explained only in 
part, if at all, by the separate theories of the strong and the electroweak interactions. Some 
of them seem to invite a further unification of the strong, weak and electromagnetic 
interactions [131]. Among the hints of deeper patterns is the striking resemblance of quarks 
and leptons. Particles in both groups are structureless at current experimental resolution. 
Quarks possess color charges, whereas leptons do not, but both carry a half unit of spin 
and take part in electromagnetic and weak interactions. Moreover, the electroweak theory 
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itself suggest a relatian between quarks and leptons. Unless each of the three lepton 
families (the electron and its neutrino, for example - see above) can be linked with the 
corresponding family of quarks ( the u and d quarks, in their three colors) the electroweak 
theory will be beset with mathematical inconsistencies. 

“In this regard, it is necessary to add interesting results of a non - accelerator study of 
the strong interaction in LiD crystals, The measured dependence of the strong interaction 
constant on the distance between nucleons in the deuterium nucleus [15] showed that it 
changes from zero (LiH crystals, the hydrogen nucleus does not have a strong interaction) 
to maximum value in LiD crystals (the deuterium nucleus has a strong interaction). From 
experimental value of the isotopic shift (0.103 eV) in the paper [132] we have obtained a 
residual strong coupling constant equal to 2.4680. 


What is known about the fundamental forces also points to a unification. All three can be 
described by gauge theories, which are similar in their mathematical structure. Moreover, 
the strengths of the three forces appear likely to converge at very short distances, a 
phenomenon that would be apparent only at extremely large energies Fig. 20. We have 
seen that the electromagnetic charge grows strong at short distances, whereas the strong, 
or color, charge becomes increasingly feeble. Might all the interactions become 
comparable at some gigantic energy? 


Fig. 20. Evolution of the convergence of coupling constants ( ae - 
electromagnetic; a,, - weak interaction ‘and a; - strong interaction) on extremely high 
energies, which are equivalent to very small scales of distance («= 10-9 cm) is 
expected in unified theories (scheme). 
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If the interactions are fundamentally the same, the distinction between quarks, which 
respond to the strong force, and leptons, which do not, begins to dissolve. In the simplest 
example of a unified theory, put forward by authors of the paper [133], each matched set of 
quarks and leptons gives rise to an extended family containing all the various states of 
charge and spin of each of the particles. In mathematical consistency of the proposed 
organization of matter is impressive. Moreover, regularities in the scheme require that 
electric charge be apportioned among elementary particles in multiples of exactly 1/3, 
thereby accounting for the electrical neutrality of stable matter. The atom is neutral only 
because when quarks are grouped in threes, as they are in the nucleus, their individual 
charges combine to give a charge that is a precise integer, equal and opposite to the 
charge of an integral number of electrons. If quarks were unrelated to leptons, the precise 
relation of their electric charges could only be a remarkable coincidence. 

In such a unification only one gauge theory is required to describe all the interactions of 
matter. In a gauge theory each particle in a set can be transformed into other particle. 
Transformation of quarks into other quarks and of leptons into other leptons, mediated by 
gluons and intermediate bosons, are familiar. A unified theory suggests that quarks can 
change into leptons and vice versa. As in any gauge theory, such an interaction would be 
mediated by a force particle: a postulated X or Y boson. Like other gauge theories, the 
unified theory describes the variation over distance of interaction strengths. According to the 
simplest of the unified theories, the separate strong and electroweak interactions converge 
and become a single interaction at a distance of ~ 10° centimeter, corresponding to an 
energy of ~ 3 - 10'® Gev (see Fig. 20). Such an energy is far higher than may ever be 
attained in an accelerates, but certain consequences of unification might be apparent even 
in the low world we inhabit (see, also [134]). The supposition that transformations can cross 
the boundary between quarks and leptons implies that matter, much of whose mass 
consists of quarks, can decay. If, for example, the two u quarks in a proton were to 
approach each other closer than 10°29 centimeter, they might combine to form an X boson, 
which would disintegrate into positron and ad quarks. The antiquark would then combine 
with the one remaining quark of the proton, a d quark, to form a neutral pion, which itself 
would quickly decay into two photons. In the course of the process much of the proton’s 
mass would be converted into energy. 

The observation of proton decay would lend considerable support to a unified theory. It 
would also have interesting cosmological consequences. The Universe contains far more 
matter than it does antimatter. Since matter and antimatter are equivalent in almost every 
respect, it is appealing to speculate that the Universe was formed with equal amounts of 
both. If the number of baryons - three - quark particles as the proton and the neutron, which 
constitute the bulk of ordinary matter - can change, as the decay of the proton would imply, 
then the current excess of matter need tot represent the initial state of the Universe. 
Originally matter and antimatter may indeed have been in equal quantities, but during the 
first instants after the big bang, while the Universe remained in a state of extremely high 
energy, processes that alter baryon number may have upset the balance. A number 


-419 - 


experiments have been mounted to search for proton decay. The large unification energy 
implies that the mean lifetime of the proton must be extraordinary long - 10°° years or more. 
To have a reasonable chance of observing a single decay it is necessary to monitor an 
extremely large number of protons; a key feature of proton - decay experiments has 
therefore been large scale. The most ambitious experiment mounted to date is an 
instrumented tank of purified water 21 meters on a side in the Morton salt mine near 
Cleveland. During almost three years of monitoring none of te water’s more than 
10°%protons has been observed to decay, suggesting that the proton’s lifetime is even 
longer than the simplest unified theory predicts. Income rival theories, however, the lifetime 
of the proton is considerably longer, and there are other theories in which proton’s decay in 
ways that would be difficult to detect in existing experiments. Furthermore, results from 
other experiments hint that proton can indeed decay. 

Concluding this part we note that although the SM free of inconsistence, it is incomplete; 
one left hungry for further explanation [135 - 137]. The model does not account for the 
pattern of quark and lepton masses [138] or for the fact that although weak transitions 
usually observe family lines, they occasionally cross them. The family pattern itself remains 
to be explained. Why should there be three sets of quarks and leptons? 

The history of elementary particle physics lends support to the idea that properties not 
calculable in the SM, such as mass spectrum of quarks and leptons, can be understood by 
going to a deeper level of structure in which many of particles of the SM are composite [6]. 
According Greenberg [139] the next candidates for compositeness are the dirons. The 
relatively large number (six) of flavors of quarks and leptons and their easy arrangement 
into generations is a superficial hint that they may be composite. If dirons, that is quarks 
and leptons, are composed of preons, there is an energy (denote by A) around which dirons 
are dissociated into preons; that is A/c? is the mass scale of compositeness. The radius of 
a composite diron is of order fc/A. If the currently known dirons are the ground states of a 
composite system, we can expect many excited partners of the dirons at energies of the 
order A. According the results of the paper [139] A is greater than about 1 TeV. 

As we know in the SM, the quark and lepton masses are proportional to the quark and 
lepton Higgs couplings. This origin of mass Is likely to carry over to preon models. Then 
quark and lepton masses would be associated with the composite vertex functions 
connecting the Higgs particles and the quarks and leptons. Previous experience has been 
that mass spectra are calculable in terms of a deeper level of structure: The Bohr model of 
atom, the shell model of nuclei and the quark model of hadrons have all been highly 
successful in this regard. These examples lend strong support to the hope that a preon 
model of quarks and leptons will determine the mass spectrum of quarks and leptons (the 
more details see [140)}). 


4. Experimental manifestation of Quantum Chromodynamics. 


All the below non - acceleratory experimental results were obtained on compounds 
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distinguished one neutron in the nucleus. Here, the new degree of freedom is the mass of 
the isotope. LiH and LiD crystals have the largest mass isotope effect. Natural diamond #C 
and synthetic #C also differ in a single neutron. This also applies to graphene. 

Isotopic substitution only affects the wavefunction of phonons; therefore, the energy 
values of electron levels in the Schrédinger equation ought to have remained the same. 
This, however, is not so, since isotopic substitution modifies not only the phonon spectrum, 
but also the constant of electron-phonon interaction (see above). It is for this reason that the 
energy values of purely electron transition in molecules of hydride and deuteride are found 
to be different [141]. This effect is even more prominent when we are dealing with a solid 
[142]. Intercomparison of absorption spectra for thin films of LiH and LID at room 
temperature revealed that the longwave maximum (as we know now, the exciton peak [143] 
moves 64,5 meV towards the shorter wavelengths when H is replaced with D. For obvious 
reasons this fundamental result could not then receive consistent and comprehensive 
interpretation, which does not be little its importance even today. As will be shown below, 
this effect becomes even more pronounced at low temperatures (see, also [144]). 

The mirror reflection spectra of mixed and pure LiD crystals cleaved in superfluid liquid helium are 
presented in Fig. 21. For comparison, on the same diagram we have also plotted the 
reflection spectrum of LiH crystals with clean surface. All spectra have been measured with 
the same apparatus under the same conditions. 


B15; 2.05 E,eV 


Fig. 21. Mirror reflection spectra in indicated crystals at 2 K cleaved in superfluid 
helium. Light source without crystals, curve is lamp. 


As the deuterium concentration increases, the long-wave maximum broadens and shifts 
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towards the shorter wavelengths. As can clearly be seen in Fig. 21, all spectra exhibit a 
similar long-wave structure, This circumstance allows us to attribute this structure to the 
excitation of the ground (ts) and the first excited (2s) exciton states. The energy values of 
exciton maxima for pure and mixed crystals at 2 K are presented in Table 5. The binding 
energies of excitons Ep, calculated by the hydrogen-like formula, and the energies of 
interband transitions Eg are also given in Table 5. 


Table 5. Values of the energy of maxima (in meV) in exciton reflection spectra of 
pure ancl mixed crystals at 2K, and energies of exciton binding Ey, band-to-band 
transitions E, (after [43]). 


Energy,meV  LIH LiHoeeDois LiHo40Doe0 LID ®LIH (78K) 
- CBte 4950 4967. 5003 5048 4939 | 
Eos 4982 5001 5039 5082 4970 

| Bb 42 45 48 52 44 
Ey 4992—~=«<«dSO1Ds—~*S “5051 +5095 4980 | 


The ionization energy, found from the temperature quenching of the peak of reflection 
spectrum of the 2s state in LiD is 12 meV. This value agrees fairly well with the value of 
AE 2s calculated by the hydrogen-like formula. Moreover, Ey = 52 meV for LiD agrees well 
with the energy of activation for thermal quenching of free-exciton luminescence in these 
crystals [43]. 

Going back to Fig. 21, it is hard ‘to miss the growth of Ai2, [145], which in the 
hydrogen-like model causes an increase of the exciton Rydberg with the replacement of 
isotopes (see below Fig. 22). When hydrogen is completely replaced with deuterium, the 
exciton Rydberg (in the Wannier-Mott model) increases by 20% from 40 to 50 meV, 
whereas E, exhibits a 2% increase, and at 2 + 4.2 K is AEg = 103 meV. 
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Fig. 22. Binding energy of Wannier ~ Mott excitons as function of reduced mass of 
ions based on values of reduced mass of ions for °LiH, ’LiH, °LiD, 7LiD and LiT (after 
[145]). 

This quantity depends on the temperature, and at room temperature is 73 meV, which 
agrees well enough with AE, = 64.5 meV as found in the paper of Kapustinsky et al. The 
single-mode nature of exciton reflection spectra of mixed crystals LiH,Di.. agrees 
qualitatively with the results obtained with the virtual crystal model (see e.g. Elliott et al. 
[146]; Onodera and Toyozawa [147]), being at the same time its extreme realization, since 
the difference between ionization potentials (AC) for this compound is zero. According to the 
virtual crystal model, A¢ = 0 implies that AEg = 0, which is in contradiction with the 
experimental results for LIH,D-, crystals. AS was shown above [43, 69} the change in Eg 
caused by isotopic substitution has been observed for many broad-gap and narrow-gap 
semiconductor compounds. 

All of these results are documented in Table 5, where the variation of Eg, En, are shown 
at the isotope effect. We should highlighted here that the most prominent isotope effect is 
observed in LiH crystals, where the dependence of Ey = f (Cy) is also observed and 
investigated. To end this section, let us note that Ey decreases by 97 cm‘! when ’Li is 
replaced with ®Li (see Table 21 in [69]). {t was in the original work (Plekhanov et al. [143)]) 
where the exciton binding energy Ey was found to depend on the isotopic composition that 
this change in Ey was attributed to the exciton-phonon interaction (originally with LO 
phonons). The preferential interaction of excitons with LO phonons in LiH (LID) crystals was 
later repeatedly demonstrated in the luminescence spectra [148] and resonant Raman light 
scattering [149], which consist of a phononless line (in the former case) and its LO 
repetitions. 

At present time we know that the main consequences of interaction of electron and hole 
interaction in excitons with polarization vibrations are the static screening of the lattice 
charges (introducing ¢9) and the change in the effective masses of the particles. Both these 
effects of electron-(hole)- phonon interaction can easily be taken into account, and lead to a 
change in the exciton Rydberg Ey. These corrections do not destroy the hydrogen-like 
structure of the exciton spectrum, At the same time, the non-coulombic corrections to the 
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electron-hole Hamiltonian modify the hydrogen-like structure -~ removing, for example, 
degeneration of levels with respect to orbital and magnetic quantum numbers (see e.g. 
[150] and references therein). The very fact, however, that the problem of renormalization of 
energy spectra of Wannier - Mott excitons does not admit an exact solution even in the 
limiting cases, often gives rise to a situation in which there is no agreement between the 
results obtained by different authors. Starting with the classical works of Haken [34], all 
papers may be divided into two broad classes depending on how they take cleal with the 
Coulomb interaction: between ‘bare’ electrons and holes, or between electrons and holes in 
the polaron state. In other words, first the interaction of band electrons and holes with LO 
optical phonons is taken into account, and then of the Coulomb interaction between 
electrons and holes clad in the ‘polarization coats’ is considered. As will be shown below, 
the study of exciton-phonon interaction in crystals with isotopic effect not only provides 
entirely new information, but also allows us to reconstruct experimentally the values of 
Frohlich and Coulomb interaction constants. From Fig. 22 we see that when hydrogen is 
completely replaced with deuterium, the binding energy of the exciton exhibits a 20% 
increase from 42 to 52 meV [144]. It is easy to see that in the model of virtual, crystal the 
binding energy of the exciton in LIT crystals [69] must be equal to 57 meV (see Fig. 22). 
Hence it follows that in the linear approximation the isotopic dependence of binding energy 
of Wannier - Mott excitons may be expressed as 

Ep =En(0)(1+7). (98) 

where Ep(0) is the purely coulombic binding energy of the exciton (in the frozen lattice), 
which in our case is equal to 31.5 meV, and the angular coefficient is B = 12.18 meV/M, 
where M is the reduced mass of ions of lithium and hydrogen (deuterium, tritium) ions; y = 
BM/Ev(O (see also [69]). From the standard equation for the Coulomb binding energy of 
the exciton 


Ep= <4 (99) 


"pee, : 
we get the dimensionless constant of Coulomb interaction: 
qn? = 482 = 0.47. (100) 
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Comparing the value of ?= 0.47 and the constant of Fréhlich exciton-phonon interaction 
g? =0.33 [59] we see that they are close enough. This implies that both the Fréhlich and the 
Coulomb interactions between electrons (holes) and LO phonons in exciton must be treated 
with equal attention. 

Under the continuous optical excitation a stationary population in the exciton states can 
be created. Due to the free motion in the crystal and the interaction with the crystal lattice 
the gas of free excitons is spread over some region of kinetic energy [86]. However, usually 
only the lowest n = 1S (ground state) exciton is populated (Fig. 4) at low temperatures [151]. 
The photoluminescence emission of free excitons can take place either in resonance with 
the exciton absorption line (a so - called zero-phonon luminescence) or can be shifted in 
energy due to the simultaneous creation of phonon (phonon - assisted luminescence). In 
the resonant excitons emission only excitons with small wavevectors of the order of the 


photon wavevector K ~0 can participate. Phonon - assisted luminescence of free excitons 
in polar compounds is mainly due to the creation of longitudinal optical (LO) phonons (151, 
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36]. Since the energy of LO phonons has a weak dependence on wavevector @ (see, e.g. 
[152], the spectrum of emitted photons is simply related to the initial distribution of exciton 
energy. Fig. 23 shows the free excitons photoluminescence in LiH and LiD crystals at 2 K 
with photoexcitation above the intrinsic absorption edge [151]. 


Fig. 23. Free exciton emission spectra at 2 K in LiH and LiD crystals cleaved in 
superfluid helium. 


The free exciton photoluminescence spectrum of LiH crystals cleaved in superfluid liquid 
helium consists of the narrow zero - phonon line and its wider phonon replicas 
corresponding to the radiative annihilation of excitons accompanied by excitation from one 
to five LO phonons. The zero - phonon emission line is almost in resonance with the 
reflection line in the exciton ground state (E,,_ 3 = 4.950 eV for LiH crystals at 2 K [151, 
43]), which is due of the direct electron transition Xi - X4 of the first Brillouin zone [52]. 
Phonons replicas form an equidistant series to the red from resonance exciton emission 
line. The difference in energies of these replicas, as early [145] is 140 meV which is close to 
the energy of LO phonons in the center of the Brillouin zone obtained in [154, 155]. The 
photoluminescence spectrum of a LiD crystals with a pure surface shown also in Fig. 23 is 
similar in many respects to the free excitons luminescence spectrum of LiH crystals. 
However, there are some differences. These differences are related with the next 
manifestations: 

1. The short - wavelength shift as whole of the free excitons photoluminescence 
spectrum of LiD crystals on 103 meV relatively the spectrum of LiH crystals. 

2. In the case of LiD crystals the energy difference between lines in the spectrum is on 
average 104 meV corresponding to the energy of the LO phonons in the I - point of the 
Brillouin zone [43, 151, 153, 154). 

Firstly the zero-phonon emission line of free excitons in LiD crystals shifts to the short - 
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wavelength side on 103 meV. These results directly show the violation of the strong 
conclusion (see, e.g. [5, 16]) that the strong force does not act on leptons. The second 
difference concludes in less value of the LO phonon energy, which is equal to 104 meV. 
When light is excited by photons in a region of fundamental absorption in mixed LiH,D:-. 
crystals at low temperature, line luminescence is observed (Fig. 24), like in the pure LiH and 
LiD crystals. As before [145], the luminescence spectrum of crystals cleaved in superfluid 
liquid hellum consists of the relatively zero - phonon line and its wide LO replicas. For the 
sake of convenience, and without scarfing generality, Fig. 24 shows the lines of two 
replicas. Usually up to five LO repetitions are observed in the luminescence spectrum as 
described in detail in [43]. In Fig. 24 we see immediately that the structure of all three 
spectra is the same. The difference is in the distance between the observed lines, as well 
as in the energy at which the luminescence spectrum begins, and in the half - width of the 
lines. 


5.00 480 E.eV 


Fig. 24. Photoluminescence spectra of free excitons in LiH (1), LiH,D; (2) and 
LiD (3) crystals cleaved in superfluid helium at 2 K. Spectrometer resolution is 
shown. 


The simplest approximation, in which crystals of mixed isotopic composition are treated 
as crystals of identical atoms having the average isotopic mass, is referred to as virtual 
crystal approximation (VCA) [156]. Going beyond the VCA, in isotopically mixed crystals 
one would expect local fluctuations in local isotopic composition within some effective 
volume, such as that an exciton. As follows from Fig. 21, excltons in LiH,Dj-, crystals 
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display a unimodal character, which facilitates the interpretation of their concentration 
dependence. As will be shown in the discussion the concentration dependence of the power 
of strong nuclear interaction, i.e. dependence on the neutron concentration. 

Because of the low intensity of scattered light, and thanks to the high resolution of 
modern spectroscopic instruments, the development of high sensitive techniques of 
detection of weak optical signals (photon counting mode, optical multichannel analyzers, 
optical linear arrays and the specialized systems [157]), the light scattering method has 
become one of the most common technique for the research. 

In the last five decades, inelastic light scattering spectroscopy has become one of the 
most powerful and widely used optical method for obtaining quantitative information about 
the interaction of light with matter. For example, it has been used to investigate resonant 
Ramans cattering (RRS) in crystals. For a displacement of the excitation like frequency 
toward long wavelength by an amount compared to the exciton resonance, e.g. E: ( E,-4s 
(where E,..;, Is the energy of the exciton ground state) intense light scattering is observed 
[154]. The spectrum of this scattering is shown in Fig. 25. As in the luminescence, the 
process of energy relaxation take place, mainly with emission of LO phonons. This is shown 
by the character of the structure in the scattering spectrum. Indeed, as for the case of the 
excitation well within the exciton zone (see Fig.4), the energy difference between the peaks 
in the scattered spectrum equals the energy of the LO phonon in the F - point of Brillouin 
zone [153 - 155). 
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Fig. 25. Resonant Raman scattering spectrum of a LiD crystals cleaved in liquid 
helium at the excitation E = 4.992 eV at 4.2K. The arrow is indicates the energy 
position of the ground exciton states. 
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. The relatively large half - width of the scattered peaks should be noted [149]. Additional 
investigations have shown that their half - width are always larger than that of the excitation 
line [43]. In addition to the LO(Ir) lines the RRS spectrum contains two more bands the 
maxima of which are displaced by twice the energy of the TO(I) and LO(X) phonons from 
the exciting line (for LIH foto) = 76 meV and harojx) = 117 meV [43, }]. More precise 
measurement of the dependence of the second - order Raman spectra for LiH [158] and 
LiD on the excitation energy (below E,, _. ;,) is shown in Fig. 26. The presence of the second 
- order TO(I) phonons fn the RRS spectrum may be explained by a relatively strong 
scattering deformation mechanism in these crystals, where, however the main mechanism 
as was seen from Figs, 24 and 25, is Frohlich mechanism of intraband scattering [158].The 
longwavelength displacement of the excitation line frequency relatively exciton resonance a 
monotonic decrease the intensity of RRS spectrum as whole more than 60 fold in both LiH 
and LiD crystals (see, also Fig. 26). 
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Fig. 26. The dependence of the shape and maximum frequency of the 2LO(T) 
Raman line on the excitation energy of LiD crystals. 


As was indicated above another very interesting example is carbon. Carbon, one of the 
most elements in nature, still gives a lot surprises. It is found in many different forms - 
allotropes - from zero dimensional fullerene, one dimensional carbon nanotubes, two 
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dimensional graphene and graphite, to three dimensional diamond - and the properties of 
the various carbon allotropes can vary widely {62, 63,65]. Fullerenes are essentially hollow 
carbon shells of various sizes. the most well - known of these is a 60 - carbon unit called 
buckminster fullerene or Cg (more details see below). For instance, diamond is the hardest 
material, while graphite is one of the softest; diamond is transparent to the visible part of 
spectrum, while graphite is opaque; diamond is an electrical insulator, while graphite and 
graphene are a conductors. Very important is that all these different properties originate 
from the same carbon atoms, simply with different arrangements of the atomic structure. 
Below we describe the new phenomena of the carbon - Isotope effect in diamond as well as 
graphene. Crystals #C and C diamond differ only one neutron. 

Due to the indirect gap of E, = 5.47 + 0.005 eV (295 K), at K = 0.76 X, diamond has 
intrinsic phonon - assisted free exciton luminescence lines (see, e.g. [68, 69]). The change 
of the indirect gap of diamond between pure '*C and "C crystals has been determined by 
Collins et al. [159]. The luminescence spectra of the natural ('*C) and synthetic ('°C) 
diamond at electron excitation were investigated by Collins et al. [159], Ruf et al. [160], 


Watanabe et al. [161]. Fig. 27 compares the free exciton luminescence for a natural 
diamond with that for a synthetic diamond. 
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Fig. 27. Spectra measured at 77K phonon - assisted free cathodoluminescence 
feature (A, B and C) and the phonon assisted bound - exciton feature (D) from natural 
semiconducting '*C diamond and a °C synthetic diamond (after [159}). 


The peaks labeled A, B and C are due, respectively, to the recombination of a free 
exciton with the emission of transverse - acoustic, transverse ~ optic and longitudinal - optic 
phonons having wavevector + Kmin and quanta (in '2C (43, 114]. 
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hata = 87 + 2; hato = 141 + 2; hay.o = 163 + 1 meV. (101) 


Features B, and Bs are further free exciton processes involving the above TO phonon 
with one and two zone - center optic phonons respectively. As we can see from Fig. 27 the 
isotope shift of the free exciton luminescence spectrum of °C diamond is equal 16.5 + 2.5 
meV [69]. The more detailed and quantitative investigation of E, ~ f(x), where x is the 
isotope concentration was done by Ruf et all. [160] (Fig. 28) and Watanabe et all. [161] (Fig. 
29) where five samples of diamond with different concentrations x were studied. 
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Fig. 12. Cathodoluminescence spectra of isotopically modified diamond at 36 K. 
Intrinsic phonon - assisted recombination peaks are labeled in the top spectrum, 
those from boron - bound excitons in that at the bottom. The spectra are normalized 
to the intensity of the B peak and vertically offset for clarity (after [160}). 


Watanabe et al. have concluded that the maximum change of the band gap due to 
substitution of '*C by °C Is AE, = 15.4 meV. 
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Fig. 29. Luminescence spectra of free excitons in homoepltaxlal diamond films 
grown from mixture of methane in hydrogen by means of a microwave plasma - 
assisted CVD. The spectra illustrate the effects of isotope composition '*C}3,C, (x = 
0.001; 0.247; 0.494; 0.740; ancf 0.987) mixed in the CVD gas phase. All spectra are 
normalized to the same height (after [161]) . 


This value is in good agreement with the estimate of 16.5 + 2.5 meV by Collins et all. 
[159] using a zero - point renormalization obtained from a fit of the experimental 
temperature dependence of the band gap (see, also [68)). 


It is clear that the lattice - dynamic properties of a crystal are directly affected by the 
atomic mass. To a first approximation, phonon behave like harmonic oscillators with 
frequencies [162] 


oom, (102) 


Crystals containing various isotope sare usually described within virtual crystal 
approximation (VCA) [162] where m in Eq. (102) is replaced by the average atomic mass 


i = Sccimi. (103) 
{ 


It is obtained from the sum over the isotopes masses mj; and concentration cj. In spite of 
VCA simplicity, the model describes the general feature of the lattice dynamics of mixed 
alkali - halide crystal sufficiently well (for details see [162]). Isotopes are ideal for lattice - 
dynamitic investigations of crystals. One can make use of the unique isotopic properties by 
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varying the isotopic composition. Isotope substitution helps to disentangle the individual 
contributions of anharmonic and disorder - induced effects and to clarity the origin of 
phonon broadening mechanisms. Due to the fact that substitution of the isotopic mass in 
semiconducting crystals a small variation of E, [81], perturbation theory is applicable 
(excluding LiH crystals [43)]), 

Raman spectroscopy is a powerful means to gain experimental access to phonons and 
their interaction and scattering mechanisms. All studies presented in this paragraph are 
restricted to stable, i.e. non - radioactive isotopes. About 300 stable and 1000 radioactive 
isotopes are known today. Some elements are isotopically pure (for example, Co), while 
others may contain numerous isotopic modifications (foe example, Sn has 10 stable 
isotopes with atomic masses ranging from 112 to 124, while Xe has 23 isotopes, 9 of which 
are stable (see, Table 1 in [68]) [163]. 

In this part, the modern understanding of first - order Raman slight scattering spectra in 
isotopically mixed elementary and compound (CuCl, GaN, GaAs, GaP) semicoriductors 
having a zinc blende structure is described. It is well - known that materials having a 
diamond structure (C, Si Ge, a - Sn) are characterized by the triply degenerate phonon 
states in the I" - point of the Brillouin zone (k = 0) (see, e.g. [85]). Isotope effect in light 
scattering spectra in Ge crystals was first investigated by Agekyan et al. [164]. A more 
detailed study of Raman light scattering spectra in isotopically mixed Ge crystals has been 
performed by Cardona and coworkers [70]. Watanabe et all. have concluded that the 
maximum change of the band gap due to substitution of *C by °C is AE, = 15.4 meV. 
This value per one neutron and on seven neutrons we get 15.4 x 7 ~108 meV. This value is 
very close to the observed one (103 meV) in LiD crystals. 

The cubic modifications of crystalline carbon, diamond, is characterized by a tetrahedral 
coordination underlying its structure dictated by sp? bonding between the nearest neighbor 
atoms. Diamond has two atoms per primitive (Bravais) cell. The strong covalent bonding 
and the light mass of the constituent atoms result in a large frequency (see, Eq. (102)) for 
zone center, Raman active, triply degenerate F, mode [152]. Its crystalline perfection and 
transparency [165] make diamond ideally suited for inelastic light scattering. First Raman 
study of the dependence the frequencies of phonons on the isotopic composition was 
conducted on diamond by Chrenko [167] in 1988. Several publications on diamond by other 
groups followed later [88, 90, 91, 92]. Most clear results was obtained by Hanzawa et al. 
(166] (see Fig. 30), 
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Fig. 30. First - order Raman scattering in isotopically mixed diamond crystals 


12C15C,_,.. The peaks A, B, C, D, E and F correspond to x = 0.989; 0.90; 0.60; 0.50; 0.30 
and 0.001{after [166)). 


First - order Raman light scattering spectrum in diamond crystals also includes one line 
with maximum @zto() = 1332.5 cm [170]. In Fig. 14 the first - order scattering spectrum 
in diamond crystals with different isotope concentrations is shown [166]. As was shown in 
{170], the maximum and width of the first - order scattering line in isotopically - mixed 
diamond crystals are nonlinearly dependent on the concentration of isotopes x. The 
maximum shift of this line is 52.3 cm~', corresponding to the two limiting values of x = 0 and 
X = 1, The effect of the isotopic '*C to %C ratio on the first - and second - order Raman 
scattering of light in the diamond has been investigated in [172]. As °C content is increased 
from the natural ratio (1*C/5C = (1 - x)/x, where x = 0.011 to the almost pure °C (x = 0.987) 
the whole spectrum has shifted towards longer wavelength (Fig. 15) in good agreement with 
the expected M-°5 frequency dependence on the reduced mass M. For an approximately 
equal mix of the two isotopes, the authors reported that the feature seen in the above two - 
phonon spectra were either broadened or unresolved. 
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Fig. 31. The Raman spectra of natural and a °C diamond. The spectra show the 
dominant first - order Raman - active F., line and the significantly weaker quasi - 
continuous multiphonon features (after [171]}. 


The analogous picture was obtained on the isotope dependence of Raman spectra in 
silicon [172]. Fig. 16 shows the measured frequencies of the first - order Raman active 
modes in silicon for six different isotopic composition. 
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Flg. 32. Typical low temperature Raman spectra of silicon with various isotopic 
composition (after [ 172)). 


As was shown above, in two - dimensional graphene, carbon atoms are periodically 
arranged in an infinite honeycomb lattice (see Fig. 6). Such an atomic structure is defined 
by two types of bonds within the sp? hybridization [60, 78]. Elastic and inelastic light 
scattering are powerful tools for investigating graphene [173 - 178]. Raman spectroscopy 
allows monitoring of doping, defects, disorder, chemical and isotope [68, 69] modifications, 
as well as edges and uniaxial strain. All sp? - bonded carbons show common features in 
their Raman spectra, the so - called G and D peaks (see, e.g. Fig. 8 in [179]), around 1580 
and 1360 cm~! (see, e.g. [174, 175]). The G peak (see, also below Fig. 34) corresponds to 
the E2, phonon at the Brillouin zone center.(F - point). The D peak is due to the breathing 
modes of six - atom rings and requires a defect for its activation. It comes from TO phonons 
around the Brillouin Zone K point and it is activated by an intravalley scattering process 
174]. The 2D peak is the second order of the D peak. This is a single peak in monolayer 
graphene, whereas it splits into four bands in bilayer graphene, reflecting the evolution of 
the band structure [7, 22]. The Raman spectrum of graphene also shows significantly less 
intensive defect - activated peaks such as the D’ peak , which lies at ~ 1620 cm~!.This is 
activated by an intravalley process i.e. connecting two points belonging to the same cone 
around K (see, Fig. 2) in [175]). The second order of the D’ peak is called 2D’ peak. Since 
2D and 2D’ peaks originate from a Raman __ scattering process where momentum 
conservation is obtained by the participation of two phonons with opposite wave vector (¢ 
and -g), they do not require the presence of defects. Thus, they are always visible in the 
Raman spectrum (see cited papers [173 - 178] and references therein). 

Graphene is one unique material which shows properties not fond in other materials. 
One of these unique features of graphene is the influence of long range strains on the 
electronic properties. The possibility of tuning the dynamics of carriers as well as phonons 
by appropriately designed strain patterns opens the way for novel applications of graphene, 
not possible with any other materials (see, e.g. [183] and references therein). At present 
time we have several reports, which have examined graphene properties under unlaxial 
deformation [178, 183, 184, 185). 

Strain can be very efficiently studied by Raman spectroscopy since this modifies the 
crystal phonon frequency, depending on the anharmonicity of the interatomic potentials of 
the atoms. Raman spectra of strained graphene show significant redshifts of 2D and G 
band (see Table 6) because of the elongated of the carbon ~ carbon bonds. 


Table 6. Red shift of the G and 2D bands in the Raman spectra in graphene 
monolayers under uniaxial tensile stress. 
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Ref. Shift of G (G* and G’) band cm™'/% Shift of 2D band cm™'/% E,, meV 
15 14.2 27.8 300 
25 5.6; 12.5 21 
27 10.8; 31.7 64 
a theory = 500 


The authors of the paper [184] have proposed that by applying uniaxial strain on 
graphene, tunable band - gap at K - point can be realized. First principle calculations 
predicted a band - gap opening of ~ 300 meV for graphene under 1% uniaxial tensile strain 
(Fig. 33). Thus, the strained graphene provides an alternative way to experimentally tune 
the band - gap of graphene, which would be more efficient and more contollable than other 
methods (see, above) that are used to open band - gap in graphene. 
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Fig. 33. The band ~ gap of strained graphene with the increase of uniaxlal tensile 
strain on graphene. The magnitude of gap is determined by the gap opening of 
density of states. The insert show the calculated density of states of unstrained and 
1% tensile strained graphene. The dash line and solid dot indicate the calculated 
bandgap of graphene under the highest strain (0.78 %) [18]. 


The method of the isotope renormalization of the energy of elementary excitations in 
solid very often used in last five decades and well described in the scientific literature (see, 
for example reviews {68, 69]). At now there is a large list of the paper devoted to 
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investigation of the isotope - mixed graphene [60, 177, 186 - 190]. Chen at al. [176] have 
reported the first experimental study of the isotope effect on the thermal properties of 
graphene. The thermal conductivity K, of isotopically pure *C (0.01 of 'C) graphene 
determined was higher than 4000 W/mK (approximately two times more than it in diamond 
[80, 181]) at the measured temperature T,, ~ 320K, and more than a factor of two higher 
than the value of K in a graphene sheets composed of a 50% ~ 50% mixture of !2C and 8C. 
Raman spectroscopy transferred to the 285 nm SiO2/Si wafer was performed under 532 nm 
laser excitation [176]. The G peak and 2D band positions in Raman spectra of graphene 
with 0.01%, 1.1% , 50% and 99.2% C - isotope are presented in Fig. 34. 
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Fig. 34. Raman spectra of graphene with different isotope concentration at room 
temperature [176]. 


Isotope shift of the G and 2D bands in the Raman spectra depicted on the Fig.35 [191]. 
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Fig. 35. Peak position of G and 2D bands in Raman spectra as a function of the 
concentration of °C [191]. 


As in the case of isotope - mixed diamond [68, 69] the Brillouin - zone - center optical - 
phonon frequency varies with the atomic mass M as w ~ M~? making the Raman shift for 
3C approximately (12/13)~'? times smaller than that for '*C. The experimental difference 
between the lowest 99.2% '%C and the highest 0.01% “C peak is ~ 64 cm~ which is 
according [176] in agreement with the theory, and attests for the high quality of isotopically 
modified graphene. By the way we should indicate that in the Raman spectra in diamond 
(with sp® - bond) analogous shift of first - order line in Raman spectrum is equal 52.3 cm~! 
[170], which is consistent with the isotope mass ratio. Substituting a light isotope (!2C, H) 
with a heavy one increases the interband transition energy in the case %C!C,.. on 14.7 
meV and LiH,.D,, on 103 meV [150]. Taking into account a more soft bond (sp* - bond 
instead sp? - bond in diamond) isotope - induced band - gap opening in graphene of some 
hundreds meV (up to E, of Si) was predicted in paper [60]. Such estimation of the value of 
isotopical band - gap opening in graphene agrees with not only the results of paper [170] 
but with very small value C4, = 0.5 - 10* dyn/em?. Such small value indicates on the strong 
electron - photon interaction - main reason renormalization of electron excitation energy (for 
the details, see, e.g. [69]). Very close to isotopically renormalization of electronic excitation 
energy is the hydrogenation of graphene [60, 179]. in last mechanism there is observable 
band - gap opening in graphene. We should add that use deuteriun instead of hydrogen we 
may increase the value of E, [69]. Thus, isotope substitution will be very useful method for 
renormalization of the band - gap in graphene - future semiconducting material. Moreover, 
this method allows to control not only of the strong nuclear interaction (quantum 

* chromodynamics) but taking it into account at the renormalization of the electromagnetic 
interaction (quantum electrodynamics) [111]. Adding °C makes magnetic materials isotope 
out of graphene. 
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Discussion. 


Traditionally nuclear - electron interaction (in our case neutron - electron interaction) 
taking into account the solving Schrédinger equation using Born - Oppenheimer (adiabatic) 
approximation [192]. Below we consider the essence of the Born - Oppenheimer 
approximation of the nuclear - electron interaction on the molecule example.. The 
Schrédinger equation for the molecule 


Aw(?, R) = Ev(?, R) (104) 


here ‘Y(P, R) ) is an eigenfunction and E the corresponding eigenvalue, A the Hamiltonian 
operator of molecule and 7 is electron and Ris nuclear distance. The simplest approximate 
method for solving the Schrédinger equation (104) uses, as indicated above, the so - called 
Born - Oppenheime approximation [192]. This approximation is used that the nuclei are 
much heavier than an electron and, consequently, move very slowly in comparison with the 
electronic motion. 

lf the nucleus to be fixed in some configuration for this particular R (distance between 
nuclei) we get the next differential equation 


es my Ar¥(?, Fy - & x fe.y(P, R) + V(P, RyY(?, Fl) = EY(?, FR), (105) 


where A is Laplasian [100], m and M electron and nucleus mass, index i for electrons 
and @ is index for nuclei. 
Further we put 


w(P, R) = X17, ROR) (108) 


here oR) is the nuclear wavefunction, which depends on nuclear distance R and on 
electronic state. 
After substitution (106) in (105) we have 


HOM DAN, R) - £ ps fev(?, RD(R) + WP, RYY(?, RQ) = E X(?, RM). 
(107) 
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Further we neglect the next members in (107) 


Cw(?, R)=-£> AHOR)a.K?P, R) + 2v-Q(R)V.X(?, A). (108) 


Operator C is call the non - adiabatic operator, which describe the nuclear - electron 
interaction. after the last oeration the equation (108) splits on two equations: 


“Bir DAIX(?, R) + VCP, YUP, F) = EXAX(?, Ri (109) 
and 
EL HOR) + E1RIOR) - WOR), (110) 


Where E.R) is electronic energy and W is the nuclear energy. Every equation have the 
view of the Schrédinger equation. The first equation (109) describes the electron stationary 
state in adiabatic approximation. Equation (110) is an equation for a wavefunction of the 
nuclei alone. Note, please that Equation (109) contains the nuclear charges through the 
Coulomb potential, but is doesn’t include any reference to nuclear mass. Thus E,(R), the 
eigenvalue (energy) of the electronic Schrédinger equation (109) is the same for all 
isotopes, and consequently, is isotopically independent. The independent E,(R) on isotopic 
substitution is the essence of the Born - Oppenheimer approximation. However, we must 
repeat, that the Born ~- Oppenheimer approximation [192] is the standard anzatz to the 
description of the interaction between electrons and nuclei in solids (see, e.g. [193)]). 
Moreover, comparison the non - accelerator experimental results on the luminescence 
(reflection) and light scattering [111] in the crystals which differ by a term of one neutron 
only is allowed to the next conclusions; ; 

1. At the adding one neutron (using LID crystals instead LIH ones) is involved the 
increase exciton energy on 103 meV. 

2. At the addition one neutron the energy of LO phonons is decreased on the 36 meV 
that is direct seen from luminescence and scattering spectra [58]. 

3. in the isotope effect, the energy of acoustic phonons does not depend on the 
replacement of H by D [69], which is proved by the identical structure of the light scattering 
spectra [111]. Along with this, a small change in the energy of optical (LO) phonons (36 
meV) observed in the luminescence and light Scattering spectra indicates a non - electron - 
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phonon renormalization mechanism of the energy of zero - phonon emission line of free 
excitons in LiD crystals, The same results have been jbtained for diamond. 

The X - ray diffraction investigations show that the LiH.,.Di_. mixed crystals [196] form a 
continuos row of the solid solution and testify themselves like a virtual crystal with a variable 
lattice constant (and exciton radius, respectively) that obeys Vegard’s law (see Fig. 36). 
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Fig. 36. Lattice constant in the 7Li(H,D) crystals plotted against the isotopic 
composition [76]. 


The measurements of the low temperature of reflection [111] and luminescence spectra 
of the whole series of mixed crystals is permitted to obtain the dependence of the interband 
transition (the long - range force dependence of strong nuclear interaction on the distance 
between nucleons in deuterium nucleus) energy on the deuterium concentration ( see below 
Fig. 37). Taking into account that the concentration of the isotopes is vary directly the lattice 
constant (see Fig. 36) we are plotting the dependence of the power of the force strong 
interaction on the distance between nucleons depicted on Fig. 37 This dependence has a 
nonlinear character ( compzre to Fig. 14). As can be seen from Fig. 37, VCA method (the 
straight dashed line) cannot describe observed experimental results. 
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Fig. 37. The long — range part of the force of strong nuclear interaction 
dependence on the distance between nucleons in deuterium nucleus. The straight 
dashed line is the linear dependence of the force dependence F, in the virtual model. 
The solid line corresponds to calculation using the polynom of second degree F, = 
f(r). Points derived from the reflection spectra indicated by crosses, and those from 
luminescence spectra by triangles. 


Traditionally nuclear - electron interaction (in our case neutron - electron interaction) 
taking into account the solving of Schrédinger equation using a model of Born - 
Oppenheimer (adiabatic) approximation [20, 192]. Since electrons are much faster and 
lighter than the nuclei by a factor nearly 2000 the electron charge can quickly rearrange 
itself in response to the slower motion of the nuclei, and this is the essence of the Born - 
Oppenheimer approximation. This approximation results the omission of certain small terms 
which result from the transformation(see above equation 108). As was shown above the 
eigenvalue (energy) of the electronic Schrédinger equation depends on the nuclear charges 
through the Coulomb potential, but it doesn’t include any references to nuclear mass and it 
is the same for the different isotopes [195]. The independent of the potential energy (the 
eigenvalue of the Schrédinger equation) is the essence adiabatic approximation. However, 
we must repeat, that the Born - Oppenheimer approximation [20] is the standard anzatz to 
the description of the interaction between electrons and nuclei in solids [193]. Now we 
should take into account a small contribution to isotopic shift through reduced electron mass 
B= ((MeMaucs)(Me + Myuct)) SO far aS Mazer is different for the hydrogen and deuterium 
nucleus. In this case the contribution equals AE ~ 6 meV. Contribution to isotope shift of the 
zero - phonon line in luminescence spectra of LiD crystals as well as Lamb shift and specific 
Coulomb potential approximately equal 1, and 1 meV, respectively. This values is less then 
isotope shift in our experiments on the two order and more. These all last results force us to 
search for new models and mechanisms of nuclear (neutron) - electron interaction including 
the results of subatomic physics, e.g. hadron - lepton interaction. 
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Qut of four known interactions, three are described by SM - the electromagnetic, weak 
and strong ones. The first two of them have a common electroweak gauge interaction 
behind them. The symmetry of this interaction SU(2), X U(1)y manifests itself at energies 
higher than ~ 200 GeV. At lower energies, this symmetry is broken down to U(1)em # U(1)y 
(the electroweak symmetry breaking): in SM this breaking is related to the vacuum 
expectation value of a scalar field [12]. The strong interaction in SM is described by the 
QCD, a theory with the gauge group SU(3)c. The effective coupling constant of this theory 
grows when the energy decreased {16]. As a result, particles which feel this interaction 
cannot exist as free states and appear only in the form of bound states called hadrons (see, 
also [5]). Most of modern methods of quantum field theory work at small values of coupling 
constant,as, [127 - 130], that is, for QCD, at high energies. Quarks and leptons, the so - 
called SM matter fields, are described by fermionic fields. Quarks take part in strong 
interactions and compose observable bound state hadrons. Both quarks and leptons 
participate in the electroweak interaction. The matter fields constitute three generation: 
particles from different generation Interact identically but have.different masses (see, e.g. 
[5]). For the case of neutrino, Yukawa [18] interactions are forbidden as well, so neutrinos 
are strictly massless in SM (see, however [194]). The gauge bosons, which are carriers of 
interactions, are massless for unbroken gauge groups U(1)e (electromagnetism - photons) 
and SU(3)c¢ (QCD - gluons), masses of W* and 2° bosons are determined by the 
mechanism of electroweak symmetry breaking. !t should noted that the forces between the 
quarks must be long - range, because the gluons have zero mass. This does not imply that 
forces between hadrons are also long range, because hadrons have zero color charges 
overall. The forces between the colorless hadrons are the residues of the forces between 
their quark constituents, and cancel when the hadrons are far apart. 

Returning to our non - accelerator experimental results, we should underline that in this 
paper we measure the strong coupling strength in crystals which differ by term of one 
neutron from each other. When we add one neutron in the hydrogen nucleus, we artificial 
activation of the strong Interaction. As far as the gravitation, electromagnetic and weak 
interactions are the same in both kind crystals (LiH and LiD), it only changes the strong 
interaction. Therefore a logical conclusion is made that the renormalization of the energy of 
electromagnetic excitations (isotopic shift equals 0.103 eV) is carried out by strong nuclear 
interaction. The short - range character of the strong interaction of nucleons does not 
possess direct mechanism of the elementary excitation (electrons, excitons, phonons) 
energy renormalization, which was observed in our low temperature experiments. Second 
reason that the interpretation of our experimental results is very difficult task because they 
are first demonstration of the violation of the strong conclusion in nuclear and particles 
physics that the strong nuclear force does not act on the colorless leptons (see, e.g. [5, 
16]). Moreover we have some contradictions taking into account that the forces between 
quarks must be long - range, because the gluons have zero mass. But as was mentioned 
above, the short range when forces between the colorless hadrons are the residues of the 
forces between their quark constituents, and cancel when the distance between hadrons is 
more than nuclear size [15]. We can see that the nuclear size transforms long - range 
interaction in the short - range strong one. It is very old question which up to present time 
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has not any theoretical explanation. 

In spite of above discussion, at present time we can distinguish the following 
mechanisms of the isotopic shift zero phonon line: 

1. Long - range electric field of the neutron’s quarks. This mechanism owing to the 
confinement quarks is limited by the boundary of the neutron. 

2. The possible new structure of the quarks and leptons - so - called preons [15, 197 - 
202). 

3. The most likely mechanism of the neutron - lepton interaction is connected to the 
magnetic - like strong field of neutron ’s quarks. Taking into account anomalous magnetic 
moment of the neuron in the paper [132] was obtained the value of strong coupling constant 
as = 2.4680 [195]. Quite large value in comparison with the accelerator technique value 
as(Mz) = 0.1198 [127]. The large value a, is thus justified to think that residual strong 
forces acting beyond nucteon could exist. a possible interpretation is to assume that in 
addition to the 8 gluons predicted by QCD SU(3). group there is a ninth gluon color singlet 
[5] 


— 4. frre a 
go= ap (it + 99 + bb) (111) 


This massless photon - like gluon may be strongly interacts between nucleons 
(neutrons) and leptons (electrons) (see, also [5, 16]. Returning to Fig. 37 we can note that 
our measurements permit to obtain value of strong coupling constant from as = 2.4680 
(pure LID crystals) to as = 0 (pure LiH crystals). Moreover, in Fig. 38 we show the 
dependence of as on neutron’s number in different substances. We can see as early in the 
case of pure LID crystals we have non - of as on the neutron’s number in different 
substances. 
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Fig. 38. The dependence of the strong force on the number of neutrons in 
different substances. 


Thus, the tentative interpretation of describing non - accelerator experimental results 
does not find consisted explanation at the change of strong interaction leaving to another 
mystery of SM. We should remind that intrinsic contradiction of Standard Model is already 
well - known. Really, the Lagrangian of QCD (theory of the strong interaction) has the next 
form (see, e.g. [35]): 


q 
where 


Vu = Ou > iQ4-AR, 
Gil, = O,AR > OnAR + gfMARAS. (113) 
‘v8 and Aj are quark and gluon fields, a=1,2,3,...8 are color indices , 2" and f"™ are Gell 


- Mann matrices and f symbols, mg - are bare (current) masses, q = u, d, §, ¢, ... different 
quarks. lt is common place [5, 16] that the Lagrangian (112) contains the members which 
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describe both free motion and interaction between quarks and gluons, which is defined by 
the strerigth couple g. Spacing of which it is necessary to remark that although the 
Lagrangian (112) possesses rather attractive peculiarities (see, also [5, 16]), its eigenstates 
are the quarks and the gluons which are not observed in free states (83, 203]. The observed 
hadrons in the experiment don't eigenstates in QCD. It is obvious to expect that the modern 
theory of QCD should finally overcome these difficulties. 

The next step our discussion is devoted to the neutron ~ electron binding energy. From 
our non - accelerator experimental results is followed that this value equals 0.106 eV on 
LIH,D;., and 0.107 eV on the diamond. In electrostatic model, it is shown [204] the neutron 
~ electron binding energy is less than analogous proton with an electron (13,6 eV) by 
approximately 128 times. Hence it follows that the theoretical value of the binding energy of 
a neutron with an electron is 106, 7 meV. Our experimental value of a neutron - electron 
binding energy is in excellent agreement with theoretical value. 

The remarkable properties of the graphene promise unheard of applications, especially 
in the semiconductor state. Especially graphene nanostructures are promising candidates 
for future nanoelectronics and solid - state quantum information technology. As was shown 
above, in the vicinity of K - points (Fig. 8), the low - energy electron/hole dispersion relation 
is proportional to momentum, rather than its square (see, e.g. [80, 205]}). This is analogous 
to the energy dispersion relation of massless relativistic electrons, so the electrons/holes of 
graphene are described as Dirac-fermions having no mass. In a simple neighbor model 
graphene is a semimetal with zero - overlap between valence and conduction bands. In 
order to make graphene a real technology, a special issue must be solved: creating an 
energy gap at K - points in the Brillouin zone (see, Fig. 5). 
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Fig. 39. The energy dispersion of semimetal (solid line) and semiconductor 
(dashed line). 


Different attempts have been made by researches, such as patterning graphene into 
nanoribbon, forming graphene quantum dots, making use of multilayer graphene sheet and 
applying an external electric field (see, e.g [80, 205]. It was shown that the uniaxial strain 
can open a band - gap in a metallic carbon nanotubes as well as carbon nanoribbon . As an 
example Fig. 40 shows the dependence of the energy of interband transition on the width of 
the graphene nanoribbons. The observed dependence is well described by the following 
expression E, = a*/ (W - W*), where a* = 0.38 eV/nm, and W* = 16 nm accounts for 
inactive edges [206]. This formula gives value of E, equal to several eV, if W is reduced to 
a size , for example, !°C, but it is impossible. The isotopic transformation of graphene, a 
semimetal into emiconductor deserves special attention [207, 208] This means of the 
cones dispersion pattern E(k) at the K point of the Brillouin zone of graphene reformes into 
the usual parabolic dispersion of band structure pattern (see Fig. 39). 
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Fig. 40. The dependence of E, as a function of width of the graphene 
nanoribbons (after [206]). 
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According to the results of Ref. 207 this method can be used to obtain E, of several 
hundred meV. Since value of E, depends on the concentration of %C in graphene, this 
mechanism is easily realized either chemically or technologically. With the opening E, of 
massless fermions (leptons = electrons) acquire mass. Estimates of different ways of 
opening E, in graphene give the following equality: E, = 2m [208]. When E, = 1 eV [207] 
leptons acquire a mass equal 0.5 eV - the value is not large, but finite. One of the 
mechanisms of acquire mass in elementary particle physics is considered in the paper 
[209]. This paper predicts that both quarks and gluons acquire running mass distribution in 
QCD, which are large at infrared momenta. The current quark of perturbative QCD evolves 
into a constituent quark as its momentum becomes smaller. The constituent quark mass 
arises from a cloud of low - momentum gluons attaching then selves to the current quark. 
This is DCSB: an essentially nonperturbative effect that generates a quark mass from 
nothing; namely it occurs even in the chiral limit. DCSB namely the generation of mass from 
nothing is a theoretically established nonperturbative feature of QCD. The solution of a gap 
equation shows that gluons are cannibals: they are particle species whose members 
become massive by eating each other (the details see [207 - 210). 


5. Conclusion. 


The artificial activation of the strong nuclear interaction by adding one (two or more) 
neutrons in atomic nuclei leads it to the direct observation of the strong interaction in low - 
temperature optical spectra of solids. This conclusion opens new avenue in the investigation 
of the constant of strong nuclear interaction in the wide value range by means the 
condensed matter alike traditional methods. The observation of an isotopic shift (0.103 eV) 
of the zero-phonon line of free excitons in the luminescence spectra of LiH (without strong 
interaction in the hydrogen nucleus) and LiD (with strong interaction in the deuterium 
nucleus) crystals was the first and direct evidence of the long-range interaction of the 
Yukawa potential. Indeed, in both crystals, the lithium ions, the proton and the electron are 
the same and, therefore, the gravitational, electromagnetic and weak interactions are the 
same, and the addition of a neutron, according to Yukawa, a strong interaction appears, the 
influence of which manifests itself in the isotopic shift. These experimental results 
demonstrate the neutron-electron binding energy (0.106 eV) which is in excellent agreement 
with the theoretical Breit estimate of 0.1067 eV Experimental observation of the 
renormalization of the elementary excitation energy of solids by the strong nuclear 
interaction stimulates its count in the process of description of the elementary excitations 
dynamics in quantum electrodynamics. Present review continuous to develop between 
nuclear and condensed matter physics. 
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Abstract 


The aim of this study is to study the Differential Cross Sections (DCS) in the presence of a linear 
polarized laser field for Bessel functions of various orders, with oxygen molecules as the target 
within a Lennard-Jones potential framework. The developed model is theoretical in which different 
assumptions are used, born approximation, Bessel function, central approximation, etc. the 
developed model ‘was computed for analysis using online MATLAB. The observation shows 
projected electron gain and loss energies in the presence of the oxygen molecules, A pronounced 
deep peak in the results suggests an alignment of incidence and final energies, potentially indicating 
new particle formation. Additionally, a minimum DCS point signifies an intense interaction and 
scattering between the incidentelectron andthe target, suggesting the formationof particles smaller 
than both the target and the election. Notably, the zero-order Bessel function exhlbits ahigher DCS 
than its higher-order counterparts, likely owing to distinct mathematical and physical properties. 
At a zeto-radian scattering angle, the DCS achieves its minimum, highlighting a heightened 
probability of particle formation due to close alignment and collision between incident and 
scattered particles. This finding underscores the significance of this pafticular angle in 
understanding particle interactions. This study's revelations offer valuable implications for 
understanding energy exchange and particle formation processes in similar contexts. Further 
research in this area could delve into specific mechanisms underlying the observed phenomena, 
potentially leading to advancements in particle interaction dynamics. Moreover, these findings 
contribute to the broader understanding of the fundamental physics governing such interactions, 
with potential applications in various fields of science and engineering. 


Kseywords: Differential cross section, Bessel function, oxygen molecules, Lennard-Jones 
Potential. scattering angle 
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IL. Introduction 


The phenomenon of evaporation and condensation at the liquid/vapor 
interface is common yet poorly understood. To delve into these processes at 
the molecular level, extensive molecular dynamics simulations have been 
conducted on Lennard-Jones (LJ) fluids, comprising monomers, dimers, or 
wimers. These simulations revealed that when LJ monomers touch a vacuum, 
evaporation rates are notably high due to strong evaporative cooling and a 
density differential near the interface. However, when the structure is reduced 
to dimers, the evaporation rate significantly decreases. The study also 
involves measuring evaporation and condensation coefficients, as well as 
velocity distributions of evaporated molecules, comparing the results with a 
model based on kinetic theory of gases. The research underscores the 
importance of considering physical quantities and appropriate parameters for 
computational material production. 


The LJ equation proves valuable in predicting material characteristics, 
facilitating interactions between pairs of neutral atoms and molecules. These 
interaction potential values can be theoretically calculated through quantum 
mechanics or confirmed through classical simulation, which can further be 
used to fine-tune experimental data. Overall, this approach allows for the 
anticipation of material properties through simulations tracking atom 
movements within a molecular dynamics framework. Lennard-Jones 
Potential (LJP) studies are crucial for analyzing the differential cross-section 
of shielding materials, especially due to the impact of atomic binding on 
radiation shielding. Materials with appropriate wavelength and bond length 
are considered ideal for shielding, as they effectively absorb and protect 
against various radiation hazards. Additionally, LJP plays a significant role 
in scattering, influencing the interaction of photons and atoms. This proves 
instrumental in studying energy, momentum, differential cross-section, and 
particle formation—a valuable method for comprehensive analysis {1-3]. 


The investigation of Lennard-Jones Potential (LJP) interactions among 
various atoms and molecules using a mathematical model and MATLAB. It 
explores the Lennard-Jones Potential direct correlation function (LJPDCF) 
for carbon-carbon, carbon-oxygen, oxygen-oxygen, water-water, silver- 
silver, gold-silver, gold-carbon, and silver-carbon interactions. The research 
unveils significant findings, such as a maximum LJP of -133.75 kcal/mole 
for silver-silver at 3.5 A and a minimum potential of -2.08 kcal/mole for 
carbon-carbon at 3.8 A. These insights have implications for molecular 
physics, computational chemistry, and molecular modeling applications [4]. 
The pivotal role of interatomic potentials in material physics models, 
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particularly those derived from first principles. These potentials underpin 
diverse phenomena like phase transitions, diffusion, and adsorption. 
Parameters are deduced from thermodynamic properties and crystal 
structures through simulations, aligning macroscopic preperties with 
experimental data. Simplified pair interaction potentials, characterized by 
v(rij), offer crucial approximations. However, limitations arise from 
neglecting electron density redistribution. Despite this, they serve as vital 
starting points for more complex potential constructions. Yet, a research gap 
emerges: the evaluation of potentials in simulating real physical phenomena 
does not sufficiently consider the impact of electron density redistribution, 
warranting further exploration [5]. 


The Lennard-Jones potential quantifies the bond interaction energy between 
two atoms or molecules, contingent on their separation distance. Utilizing 
molecular dynamics simulations, this research determines material 
properties, focusing on melting points. Employing the LAMMPS program, 
the study normalizes potential parameters by comparing simulated and 
experimental melting points. The outcome is a Lennard-Jones potential 
parameter value that provides a concise description of material melting 
points, particularly for Fe, Pb, Ni, and Cr [6]. The potential is increasingly 
used to determine the intermolecular parameters for simple molecules (H2, 
02, N2,CO, CH4), in which we add the inter-action terms in some cases 
potential iodide(K]) molecule [7]. 


0: 12 ij 6 
rs“) © 


7j=Internuclear separation, oj; =equilibrium distance and ¢,j;=binding 
energy [8]. In 1897, J.J. Thompson's discovery of the electron revolutionized 
physics beyond classical principles. Rutherford's 1911 experiment, scattering 
alpha particles off foil, led to the nucleus's revelation and birthed the 
planetary atomic model. Bohr's 1912 hydrogen model, rooted in Coulomb 
force and angular momentum quantization, ushered in new physics laws. 
Compton's 1923 experiment confirmed photons' existence, supported by 
Bothe and Geiger's validation [9]. 


The study focuses on inelastic electron-H(2s) dispersion in a linearly 
polarized laser beam, characterized by high scattering energy and moderate 
field intensities. It examines angular distributions, resonance structures, and 
states transitions in both elastic and inelastic scattering. The research 
underscores the complexity of the mathematical treatment, given the 
involvement of a projectile, atomic target, and photon. Assumptions, like 
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rapid projectile electrons and moderate field intensities, simplify the analysis. 
The Gordon-Volkov wave function delineates the projectile electron-laser 
field relationship, while first-order time-dependent perturbation theories 
describe the laser field's effect on the hydrogen atom. The first Born 
approximation handles the interaction between the hydrogen atom and fast 
projectile electrons [10]. 


The impact of circularly polarized laser fields on electron scattering by 
hydrogen atoms. It scrutinizes how the laser's polarization influences 
hydrogen and subsequently affects differential scattering cross sections. 
Dispersed electrons with 100 eV kinetic energy allow for a first-order Born 
approximation treatment. Yadav et al. observed that, under fixed laser and 
electron conditions, the differential scattering cross-sectional area decreases 
with increasing scattering angle. The polarization of the laser field 
significantly influences the differential scattering cross section [11]. The 
elliptically polarized laser field to analyze electron scattering by hydrogen. It 
examines how laser field polarization influences hydrogen atoms and 
subsequently impacts differential scattering cross-sections. Utilizing the first 
Born assumption and the Volkov wave function, scattered electrons possess 
high kinetic energy (>3000 eV) under moderately intense laser beams. The 
differential scattering cross-section area expands with increased incident 
electron energy, while changes in polarizing angle have no effect. Laser 
power and incident energy in a linearly polarized field determine the 
differential scattering cross-section in elliptical polarization, particularly in 
the highly energetic zone [12]. 


Dhobi et al. investigate inelastic scattering, focusing on low laser fields 
(visible and UV). Their model predicts that the differential cross-section 
expands as the target absorbs energy. Upon energy emission, the cross- 
section initially drops, eventually reaching its maximum. Energy emission 
occurs at specific electron voltages. Moreover, the cross-section grows with 
increasing scattering angles [13]. The study of elliptically polarized laser 
field to investigate electron scattering by hydrogen. It explores how the laser's 
polarization impacts hydrogen and subsequently influences differential 
scattering cross-sections, Applying the first Born assumption and Volkov 
wave function, scattered electrons possess high kinetic energy (>3000 eV) 
under a moderately intense laser beam. The area of differential scattering 
cross-section expands with increasing incident electron energy, with 
polarizing angle changes having no effect. Differential scattering cross- 
sections in elliptical polarization at high energy levels are determined by laser 
power and incident energy in a linearly polarized field [14]. 
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Il. Research Gap 


The research presented in the provided sources primarily focuses on various 
aspects of scattering experiments involving electrons and atoms in the 
presence of laser fields. While these studies investigate differential 
scattering cross-sections, the influence of laser polarization, and energy 
emission in scattering processes, none of them directly involve the 
application or utilization of the Lennard-Jones potential (LJP). The LJP is 
not discussed or employed in the context of these scattering experiments, 
The absence of LJP in these studies highlights a research gap in the 
understanding of how the LJP, a well-established potential energy model for 
describing interactions between atoms and molecules, might be relevant or 
applicable in laser-assisted electron-atom scattering experiments, Future 
research could explore the potential role of LJP in such scattering 
phenomena, investigating whether it can provide insights or improve the 
accuracy of theoretical models used to describe these interactions. 
Additionally, research could explore whether LJP-based models have any 
implications for predicting or interpreting experimental outcomes in this 
specific field of study. 


The motivation for further research in this field stems from a notable 
research gap identified in the existing studies. While these studies explore 
complex electron-atom scattering processes in the presence of laser fields, 
they notably overlook the application of the well-established Lennard-Jones 
potential (LJP). This omission presents an intriguing opportunity for 
investigation, as the LJP has proven valuable in describing atom-molecule 
interactions. Understanding how LJP may influence or enhance the accuracy 
of models in laser-assisted electron-atom scattering experiments could lead 
to valuable insights and potentially improve our ability to predict and 
interpret experimental outcomes. Bridging this gap could contribute to a 
deeper comprehension of these intricate interactions and advance the field 
of optical research further. 


Il. Significant 


The significance of this approach lies in its potential to bridge a critical 
research gap. While existing studies delve into intricate electron-atom 
scattering processes under laser fields, they overlook the applicability of the 
established Lennard-Jones potential (LJP). This omission opens an intriguing 
~ avenue for exploration, as LJP has shown promise in elucidating atom- 
molecule interactions. Understanding how LJP may impact or refine models 
in laser-assisted electron-atom scattering experiments could yield valuable 
insights, potentially enhancing our capacity to predict and interpret 
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experimental outcomes. This endeavor holds the promise of advancing our 
comprehension of these complex interactions, thus driving forward the field 
of optical research. 


IV. Method and Materials 


According to Bransden&Joachain, cross sections are "the ratio of the flux of 
the acting on a particle with regard to the target to the number of events of 
this kind per unit time and per unit scatterer." Cross sections can be regarded 
as transition probabilities per unit time (transition rate), per unit target 
scatterer, and per unit flux of the incident particles with regard to the target 
if the incident flux is set to unity. The quantum scattering theory's transition 
rate [15] is given by 


wi A= I Ty POE, — B) (2) 


where T;;is the transition matrix or T matrix, which is assumed to be 
integrated across energy with the density of final states in equation (2), 


mk, (21)° 


2 
Gnyae [Tri[ dQ (3) 
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The scattering particle's mass and final wavenumber are given by dQ. m and 
ky, respectively. The wave functions then normalized to a delta function, 
according to equation (3). Divided by the incident flux equation in variables 
of wavenumber, equation (3) yields the differential cross section. 
Di Ak; 


Fincident = m(2m)3 or m(2m)? (4) 


On consideration of electron density unity in a volume of (2z)?, the 
differential cross section is 


da ke (m(2n)3\7 
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The quantum free-free differential cross - sectional area can be determined if 
the T matrix for free-free transitions is discovered. We examine the 


connection between the S and the T matrices in order to determine the T 
matrix. 


Spi = Ori oo i216 (Ey = E;) Tj (6) 
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Now, S matrix for free-free scattering 
l 0 B P t 
Sp = Sn ZS ir (x/(r, IVY (r, t’))ae (7) 


where X (r, t) denotes the wave function for an electron connected to an 
external electromagnetic field and ¥ (r, t) denotes the wave function for an 
electron coupled to an external electromagnetic field while also being in the 
existence of a scattering potential V (r). 


1 Pp as Eee es 
X(r,t) = Laaarow ie (r+ | Ace)at' ivt 
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Equation (8) is Volkov wave function, E is the free electron's kinetic energy. 
Keep in mind that the momentum operator's eigenvalue is p. The Volkov 
wave function's general form is given by equation (8) to continue to 
determine 


X (r,t) 
yee ate is aie r+——sin(ot) 
= One “Fs + am ar aD NW 
22 
-i - ; sin(2ut)} (9) 


This equation explain the wave function of electrons in a laser field.The Kroll 
Watson Approximation (K WA), a common model for free-free processes, in 
this section. Since the electromagnetic field is treated conventionally yet 
quantum scattering theory is employed, the model is regarded as 
semiclassical. We use a similar reasoning as that presented by [16]. Going 
back to equation (7), we apply equation (9), 


2 Fon 
Spi = 51-5 [ (X¢ (r,t) Vr] YiCr, tt" (10) 


Using the first Born approximation, let's calculate the time integral by 
supposing that the scattered potential is low [17]. 


; pte 
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To determine, we assess the matrix element and integrate. 


+20 i l = i 
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+00 e 
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When the momentum transfer formula Q = pz — p; was utilized. Lets use 


the Bessel function with expanding plane wave Jacobi-Anger Expansion. The 
scattering amplitude of the first born approximation [15], 


m -io.- 
| eee = -.| V(rje iQ id?r (13) 
Using Bessel function and solving equation (12) we obtain, 
i +00 
qf Urovoime.e)ae 
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Now, taking the Fourier transformation of equation (13), and considering first 
imigainary and real part and also applying central approximation in real and 
imaginary on solving we get 


911 3,5 
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Now, the S matrix can be expressed from above for the n‘” element. 
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Sri 
= 5p; — 1208(Ey — 
1 2h 
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Comparing equation (6) with equation (16) and the value of fj,,,from 
equation (15), and recognizing that the T-; matrix interm of atomic unit is 


. e Q°o1t 07a" 
Tri =J, (- ha (Q. a)) (sent ses00 + |] (17) 


We determine the differential cross - sectional area for LAFF scattering for n 
photons using the solution to equation (17) in equation (5) and in atomic unit: 


dof,  256k,ne?[ Q°att_  Q3a°] (- baa - 
dQ’ —ik;~——S«* 907200" 6 Jn w (18) 
Where J is Bessel function defined [18]. Taking first term of first order 


function to remove the complixcity we have, J) = 1 and hence, DCS for zero 
order Bessel function we get 


doft, 256kyme7[ Q20"! Qa] 
do ~—siw;«~Ss«S907200—s« (19) 


Similarly, for first order Bessel function and hence the DCS for first order 
Bessel function, we get, 


(<2) - _ 256k,Q?a*n7e? cos? 8 Q2att 93g)? 55 
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Similarly, for second order Bessel function and hence the DCS for second 
order Bessel function we get 


(Z),- _ 256k," 6? Q*a4 costo[ Q2att = Q3a 5] 4 
aQ, 64k;w* 907200 6 (21) 


V. Results and Discussion 


The study employs Equations (19), (20), and (21) to analyze the Differential 
Cross Sections (DCS) for Bessel functions of zero, first, and second orders. 
Computational results, generated using online MATLAB, are depicted in 
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Figure | and Figure 2. The investigation centers on oxygen molecules within 
a Lennard-Jones potential framework. The underlying computational data is 
drawn from a reference [19]. The experiment utilizes a He-Ne laser emitting 
2 electron volts (2eV) and applies an electric field strength of 10° volts per 
centimeter (V/cm). This comprehensive approach seeks to elucidate the 
interaction dynamics of oxygen molecules under these specified conditions. 
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Figure 1: DCS vs Change in momentum 
Figure 1 provides a comprehensive representation of the Differential Cross 
Section (DCS) nature in logarithmic terms, illustrating changes in 
momentum. In this context, it's crucial to note that negative values of change 
in momentum indicate that the final energy of the electron is lower than the 
incidence energy. This signifies electron energy loss due to interaction with 
the laser field in please of target. Conversely, positive values of momentum 
change suggest that the final energy of the electron is higher than the 
incidence energy, indicating an energy gain from the laser field in presence 
of target. The prominent deep peak in the graph is particularly significant. It 
indicates instances where the incidence energy matches the final energy, 
implying a potential for the formation of new particles. This phenomenon is 
crucial in understanding the transformative potential of these interactions. 
Moreover, the deep point also serves as a notable observation. It signifies a 
minimum value of the Differential Cross Section, indicating a high level of 
interaction between the incident electron and the target. This suggests that the 
incident electron experiences a significant level of scattering, potentially 
resulting in a deviation from its initial trajectory. Furthermore, the graph 
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reveals an intriguing trend: the DCS of the zero-order Bessel function 
surpasses that of the first and second-order Bessel functions. This observation 
suggests that the fundamental Bessel function mode exhibits a greater 
propensity for interaction compared to higher-order modes. This 
phenomenon could be attributed to the distinctive mathematical properties 
and physical characteristics associated with these different orders of Bessel 
functions. 
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Figure 2: DCS with scattering angle 


Figure 2 provides a detailed view of the DCS behavior in relation to scattering 
angle. Notably, at a scattering angle of zero radians, the DCS reaches its 
minimum value. This phenomenon can be attributed to a significant physical 
process occurring at this specific angle. The minimum DCS at zero radian 
scattering angle is indicative of a high likelihood of particle formation. This 
is due to the fact that at zero radians, the incident and scattered particles are 
effectively aligned, leading to a high probability of collision. This alignment 
results in a substantial change in the trajectory of the particles, leading to a 
greater chance of particle formation. Moreover, the DCS being at its lowest 
point signifies a particularly intense interaction between the incident particle 
and the target. In practical terms, this suggests that the incident particle 
. experiences a substantial degree of scattering, meaning it deviates 
significantly from its original path. This high level of interaction is indicative 
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of the strong forces at play during this specific scattering event. The observed 
minimum in DCS at zero radian scattering angle points to a heightened 
potential for particle formation, underlining the significance of this specific 
scattering event. Additionally, the low DCS value indicates an intense 
interaction, with the incident particle being in close proximity to the target 
material. 


VI. Conclusion 


The study's findings provide valuable insights into particle interactions, 
Negative momentum changes indicate electron energy loss to the laser field, 
while positive changes imply energy gain. The occurrence of a deep peak 
suggests potential new particle formation, reflecting a high level of 
interaction with the target material. This results in significant scattering and 
trajectory alteration. Additionally, at a zero-radian scattering angle, the DCS 
reaches its minimum, indicating an elevated likelihood of collision and 
subsequent particle formation. This underscores the intensity of interaction at 
this specific angle, offering crucial understanding of the underlying physical 
processes. 
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